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ABSTRACT 

This  project  is  a passive,  water-wall  design  having  movable  insulation 
and  wood  auxiliary  heat.  The  house  is  located  about  three  miles  south  of 
Bozeman,  Montana,  The  performance  of  the  house  was  monitored  from  Decem- 
ber 1978  through  eanly  March  1979*  The  ambient  temperatures  were  abnor- 
mally low  producing  27  percent  higher  heating  degree  days. 

During  the  monitoring  period,  the  sun  provided  65  percent  of  the  heat 
load  and  the  average  house  temperature  was  19°C  (66°P),  The  wood  furnace 
provided  20  percent  and  electrical  dissipation  in  appliances  provided  15 
percent  of  the  heat  requirements.  The  collector  typically  operated  at  effi' 
ciency  greater  than  60  percent.  The  solar  system  provides  over  70  percent 
of  the  heat  requirements  on  an  annual  basis  and  about  one  cord  of  pine  is 
used  for  auxiliary  heat. 

The  disadvantages  of  this  design  are  the  daily  temperature  variations, 
average  6,4°C  (l2°P)iand  the  manual  operation  of  the  insulating  curtain  and 
the  auxiliary  furnace.  The  house  must  be  occasionally  vented  during  Sep- 
tember and  October  to  prevent  overheating. 


NOTE  ON  UNITS 

The  data  is  presented  in  metric  units.  These  units  are  in  accord  with 
the  recommendations  of  the  International  Solar  Energy  Society  as  well  as 
most  current  technical  journals.  The  overall  summaries  and  niimbers  in  the 
text  have  been  given  in  dual  units.  The  heat  load  calculations  are  pre- 
sented in  conventional  units. 

Most  readers  are  familiar  with  the  Celsius  or  Centigrade  temperature 

scale  and  are  also  familiar  with  energy  units  of  kilovratt  hours.  The  only 

other  metric  unit  which  occurs  frequently  in  this  report  is  the  Joule  (j), 

the  metric  unit  of  energy  which  replaces  the  Btu  (a  Btu  is  equal  to  1,055 

Joules),  The  mega  joule  (MJ)  is  equal  to  lO^J  (100,000  Btu/hr  = IO5  Mj/hr), 

2 

Other  abbreviations  used  in  this  report  include:  m = meters;  m = 
square  meters;  1 = liters;  h = hours;  gph  = gallons  per  hour;  and  MBtu  = 
millions  of  Btu's, 
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1.0  INTRODUCTION 

The  passive  solar  house  is  located  in  Bozeman,  Montana.  This  passive 
house  was  designed  and  built  by  Charless  Fowlkes  during  1976-77  and  it  has 
been  occupied  by  him  and  his  family  since  September  1977.  The  design  is 
described  as  a ’’passive  solar  water  wall”  and  the  auxiliary  heat  is  sup- 
plied by  a wood-burning  furnace.  The  solar  heating  system  is  integrated 
into  the  design  and  construction  of  the  house. 

2.0  DESCRIPTION  OF  HOUSE 

A photograph  of  the  house  and  the  passive  solar  collector  is  shown  in 

Figure  1.  The  house  is  a two-story  wood— framed  structure  with  a usable 
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floor  area  of  about  19O  m (2100  ft  ) . The  walls  are  framed  with  5 x 15  cm 
(2x6  in)  studs  on  60  cm  ( 24  in)  centers  and  the  roof  is  spanned  by  5 ^ 

25  cm  (2  X 12  in)  purlins.  The  walls  and  the  floor  are  insulated  with  I5 
cm  (6  in)  thick  glass  fiber  batts  and  the  roof  with  25  cm  (12  in)  thick 
batts.  A continuous  vapor  barrier  of  polyethylene  was  used  inside  all  in- 
sulated surfaces. 

Fenestration,  excluding  the  collector  glazing,  is  12  percent  of  the 
floor  area  and  all  windows  are  double  glazed.  Electricity  is  used  for 
lighting,  hot  water  and  all  appliances. 

3.0  SOLAR  COLLECTOR 

Radiation  enters  a vertical,  double-glazed  aperture  v;hich  faces  due 
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south.  The  total  area  of  the  aperture  is  43.6  m (474  ft  ).  The  center 
4 m (43  ft  ) is  not  backed  by  the  solar  wall  and  serves  as  a v;indow.  The 
solar  wall  consists  of  20  vertical  tubes  filled  with  water.  The  tubes  are 
filament  wovind  FRP,  40  cm  (I6  in)  in  diameter,  5.2  m (17  ft)  long  and  con- 
tain 674  kg  (1482  lb)  of  water  each.  The  total  storage  mass  of  the  wall  is 
about  13,480  kg  (30,000  lb).  The  tubes  are  painted  flat  green  on  the  por- 
tion facing  the  sun  and  are  spaced  3.5  cm  (I.5  in)  apart.  Heat  is  trans- 
ferred from  the  solar  wall  to  the  house  by  natural  convection  and  radiation. 

Design  calculations  indicated  that  a movable  insulation  system  v/ould  be 
needed  to  control  the  heat  loss  through  the  collector  glazing  at  night.  The 
insulation  system  selected  for  this  house  consists  of  two  thick  fabric  cur- 
tains which  resemble  window  shades.  Figures  1 and  2.  The  curtains  are 
mounted  on  a common  shaft  and  actuated  by  a geared  electric  motor.  The  out- 
side of  the  curtain  was  fabricated  from  6 ounce  white,  dacron  sailcloth  and 
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the  inside  from  6 ounce  nylon  pack  cloth  having  a urethane  coating.  The 
insulation  consists  of  four  layers  of  10  ounce  Polarguard,  a continuous 
polyester  fiber.  The  thickness  of  the  curtain  is  about  12  cm.  Sewn  along 
both  vertical  edges  of  the  curtain  is  a bolt  rope  which  runs  in  a continu- 
ous plastic  guide  attached  to  the  window  frame  to  minimize  heat  loss  by 
convection  around  the  edges. 

4.0  AUXILIARY  WOOD  FURNACE 

The  auxiliary  heat  is  supplied  by  a wood-burning  furnace,  Figure  3. 

The  hot  air  from  the  furnace  is  normally  directed  to  the  solar  storage  wall. 
The  furnace  always  operates  at  full  output  with  the  storage  wall  acting  to 
moderate  the  amount  of  furnace  heat  entering  the  house  itself. 

The  furnace  consists  of  a fire  box  and  a heat  exchanger.  It  was  de- 
signed by  the  author  and  fabricated  from  11  and  13  gauge  steel  sheet.  The 
fire  box  is  lined  with  fire  brick  and  is  a base  burning  design.  Preheated 
air  enters  the  combustion  chamber  through  three  adjustable  vents.  The  com- 
bustion gases  leave  the  fire  box  below  the  grate  and  pass  through  a heat 
exchanger.  A draft  induction  fan  is  used  when  starting  a new  fire.  Heat 
is  moved  from  the  furnace  to  the  solar  wall  with  an  air  loop  powered  by  a 
conventional  furnace  blower  and  thermostat.  The  only  hot  air  outlets  to  the 
house  are  at  the  base  of  the  solar  storage  wall.  The  cold  air  returns  from 
a manifold  located  on  the  other  side  of  the  solar  wall. 

If  the  house  temperature  drops  to  around  15°C  (60^F)  and  cold  and 
cloudy  weather  is  predicted,  a fire  is  built.  The  fire  box  is  filled  with 
about  75  kg  (165  lb)  of  pine  and  ignited.  The  furnace  reaches  full  output 
during  the  next  half  hour  and  bums  for  about  six  hours.  A second  load  of 
fuel  is  sometimes  added  to  the  fire  box  in  extreme  weather. 

A floor  plan  of  the  house  is  shown  in  Figures  4 and  5»  Note  that  the 
shop,  pantiy  and  garage  attached  to  the  north  side  of  the  house  are  not 
heated  and  thus  do  not  enter  the  heat  load  measurements. 

5.0  HEAT  LOAD 

The  calculated  heat  load  is  shown  in  Table  1.  Note  that  the  large  ex- 
panse of  glass  has  an  important  effect  on  the  heat  load  calculations.  Be- 
cause of  this  fact,  the  heat  load  is  calculated  with  the  insulating  curtain 
open  as  well  as  closed.  The  curtain  is  closed  primarily  at  night  and  open 
during  the  day.  It  is  possible  to  calculate  a weighted  aversige  heat  load 
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by  taking  into  account  average  ambient  temperatures  when  the  curtain  is 
open  and  closed.  This  weighted  average  is  also  shown  in  Table  1. 

The  shading  diagram  taken  1.5  m (5  ft)  above  groxmd  for  the  passive 
solar  collector  is  shown  in  Figure  6.  The  morning  shading  is  due  to  a moun- 
tain range  to  the  east  and  the  afternoon  shading  is  due  to  trees  west  of 
the  house.  There  is  a minimum  of  shading  at  this  site, 

6.0  TRANSDUCER  ARRANGEMENT 

The  transducer  layout  is  shown  in  Figure  7»  A total  of  27  temperature 
probes  were  used.  The  house  temperature  was  monitored  at  two  locations 
upstairs  and  two  locations  downstairs.  This  was  necessary  due  to  the  strat- 
ification of  the  air  in  this  system.  The  ambient  temperature,  shop  temper- 
ature and  crawl  space  temperature  were  also  measured. 

Four  sets  of  three  temperature  probes  v;ere  immersed  in  the  water  in 
the  centers  of  the  heat  storage  tubes.  These  probes  were  placed  at  vary- 
ing levels  in  an  attempt  to  compensate  for  stratification  and  to  arrive  at 
an  average  storage  temperature,  Figure  8.  The  temperature  change  in  the 
tubes  is  on  the  order  of  1°C  per  hour  while  a temperature  differential  of 
5*^C  may  exist  between  the  top  and  bottom  of  a tube  at  any  given  time.  An- 
alysis of  the  data  showed  that  short-term,  average  power  calculations  were 
noticeably  affected  by  the  location  and  number  of  temperature  samples. 

Daily  average  power  calculations  were  relatively  independent  of  the  probe 
locations. 

The  solar  transducer  was  positioned  behind  the  glazing  in  the  center 
of  the  solar  wall  to  measure  actual  solar  radiation  input  to  the  storage 
wall,  Figure  7-  Electric  power  input  was  monitored  by  measuring  the  cur- 
rent in  the  input  mains  with  a clamp-on  ammeter.  This  transducer  is  cali- 
brated on-site  using  a watt  meter  and  has  an  accuracy  of  at  least  3 percent. 
The  status  inputs  used  a switch  or  relay  to  record  the  status  of  the  fur- 
nace fan  and  the  insulating  curtain. 

The  furnace  output  was  determined  from  the  temperature  difference  in 
the  supply  and  the  return  air  ducts  and  the  flow  rate.  Figure  9»  The  ac- 
curacy of  this  measurement  is  about  10  percent.  Air  leaks  located  beyond 
the  measurement  station  could  further  increase  the  uncertainty  of  the  actual 
furnace  power  input.  The  weight  of  wood  burned  in  the  auxiliary  furnace  was 
determined  to  within  5 percent. 
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7.0  HOURLY  PERPORI^IANCE 

The  data  system  v;as  in  the  house  from  November  28  to  March  6.  During 
this  period  2,020  hours  of  average  data  were  recorded  for  19  individual 
channels.  The  total  data  lost  in  this  period  was  14  percent,  due  primarily 
to  one  bad  cassette  tape  in  mid-^’ebruaiy.  The  total  channel  samples  during 
the  test  were  about  25»000,000. 

The  hourly  performance  was  summarized  using  the  heat  balance  equation 
shown  in  Figure  10.  The  terms  on  the  left  hand  side  of  the  equation  were 
evaluated  using  the  hourly  data  and  formatted  by  the  computer  into  hourly 
and  daily  summaries. 

The  solar  input  consisted  of  two  terms:  the  solar  radiation  striking 
the  tubes  when  the  insulating  curtain  was  up  and  the  solar  radiation  en- 
tering the  remainder  of  the  glazing  as  simple  passive  gain. 

The  electric  power  dissipated  into  the  house  was  estimated  as  an  aver- 
age of  0.6  times  the  total  electric  consumption.  This  factor  is  less  than 
one  due  primarily  to  drain  losses  in  the  electrically-heated  domestic  hot 
water  and  an  external  vent  on  the  electric  clothes  dryer. 

The  mass  of  the  storage  wall  was  known  to  within  5 percent  and  its 
temperature  was  measured  by  the  internal  probes.  The  thermal  mass  of  the 
house  was  estimated  by  taking  into  account  the  sheetrock,  some  of  the  struc- 
tural wood  and  the  contents.  The  temperature  of  this  mass  was  assumed  to 
be  equal  to  the  average  temperature  of  the  house  air. 

The  hourly  values  of  stored  heat  changes  would  have  been  more  accurate 
if  there  had  been  a larger  sample  of  tube  temperatures.  The  hourly  change 
of  temperature  is  small  compared  to  vaxiations  existing  within  a single 
tube  as  well  as  variations  between  tubes.  This  shov;s  up  as  fluctuations 
and  phase  lag  in  the  calculated,  stored  heat  values.  Analysis  of  the  data, 
however,  shov;ed  that  these  hourly  fluctuations  average  out  and  that  the 
daily  changes  seem  accurate  and  reliable.  The  equations  used  to  form  the 
hourly  summary  are  given  in  Figure  11. 

8.0  EXPLANATION  OF  HOURLY  DATA  TABLES 

Look  at  Table  2 for  January  1,  1979  as  an  example.  The  first  column 
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is  solar  insolation  in  Mj/m  . The  next  six  columns  are  the  terms  of  the 
energy  balance  equation  of  Figure  10  and  the  units  are  mega  joules  (MJ). 

The  TUBE  STORED  HEAT  and  HOUSE  STORED  HEAT  are  positive  when  they  are  ab- 
sorbing heat  and  negative  when  they  are  delivering  heat  to  the  house.  The 
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record  starts  at  midnight  and  the  temperature  of  the  storage  wall  is  drop- 
ping about  0.4°C  (0.7°F)  each  hour  and  delivering  around  20  Mj/hr  to  the 
house.  The  ambient  temperature  drops  to  -26,5°C  (-15»7°F)  at  3 AI'l.  The 
house  temperature  is  also  dropping  and  the  thermal  mass  of  the  house  is  de- 
livering a few  MJ  each  hour.  The  house  reaches  a minimum  temperature  of 
13.5°C  (56.3°F)  at  9 The  sun  strikes  the  collector  and  the  insulating 

curtain  is  opened  aroiind  10  M (see  SOLAR  TUBE  HEAT), 

The  house  temperature  increases  during  the  day  and  reaches  a maximum 
of  21.3°C  (66°F)  at  3 PHi  just  before  sunset.  The  tubes  absorb  heat  during 
the  day  until  4 PM,  then  they  start  delivering  heat  to  the  house  through  the 
evening  and  night. 

The  hourly  heat  loss  factor  is  calculated  hy  adding  together  the  inputs 
of  the  heat  balance  and  dividing  this  by  the  temperature  difference  between 
the  inside  and  outside  of  the  house.  These  values  are  shown  in  the  CALC 
LOSS  FACTOR  column.  It  is  clear  that  the  loss  factor  increases  v;hen  the  in- 
sulating curtain  is  up  during  the  day.  The  HOUSE  TEt4P  DIFF  column  is  the 
difference  in  temperature  between  the  upper  and  lower  floor  due  to  stratifi- 
cation, The  HOUSE  TUBE  DIFF  column  is  the  storage  tube  temperature  minus 
average  house  temperature.  The  ambient  temperature,  AMBNT  TEt<IP,  reached  a 
high  of  -17.6°C  (0°F)  this  day. 

The  bottom  row  of  numbers  is  the  daily  summary  and  is  the  sum  of  all 
energy  items  and  the  average  of  the  temperatures.  The  solajr  total  input  vjas 
504  + 156  or  660  MJ,  The  electric  input  was  II4  MJ  and  the  tubes  had  a 24- 
hour  net  energy  loss  of  32  HJ. 

Appendix  I contains  a sample  of  52  days  of  hourly  data,  about  half  of 
the  data  bank.  The  format  of  this  data  is  the  same  as  discussed  in  the  pre- 
vious example.  This  data  describes  the  behavior  of  the  house  during  a 
variety  of  environmental  conditions, 

9.0  WEEKLY  SUMMARY  DATA 

Some  hourly  data  for  a very  cold  but  sunny  week  in  January  is  plotted  in 
f Figure  12,  The  solid  curves  are  solar  insolation  and  the  dotted  curves  are 

ambient  temperature,  house  temperature  and  tube  temperature.  Note  that  all 
temperatures  respond  to  the  sun.  The  solar  house,  however,  maintains  an 
ecfuilibrium  temperature  some  45°C  (80°F)  above  the  ambient  temperature  early 
in  the  week.  For  the  first  five  days  the  house  is  heated  entirely  l?y  the  sun. 
On  the  sixth  day  there  is  little  solar  radiation  and  the  auxiliary  furnace  is 
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fired,  some  of  the  heat  being  stored  in  the  tubes.  The  last  day  of  the  vjeek 
also  has  little  radiation  and  the  house  uses  stored  heat  from  the  tubes. 

Figure  13  shov;s  a bar  graph  of  the  daily  energy  use  and  distribution 
during  this  vjeek  of  January  4 to  January  10  discussed  previously.  The  top 
of  the  bar  graph  represents  the  total  heat  dissipated  by  the  house  during 
the  24—hour  period  beginning  at  midnight.  The  top  bar  element  represents 
electrical  input  and  the  shaded  area  represents  total  solar  input.  The 
last  element  represents  total  stored  energy.  If  this  element  lies  belov; 
sero,  there  v;as  a net  gain  in  stored  energy.  Average  house  and  ambient  air 
temperatures  are  plotted;  the  bars  sho;v  the  daily  extreme  values. 

During  the  first  five  days  of  this  very  cold  week,  the  sun  was  meeting 
the  heat  load  but  the  net  stored  energy  was  relatively  small.  Aiixiliary 
v7ood  heat  vjas  required  on  the  sixth  day,  v;hich  tvas  overcast.  The  following 
day  v;as  also  overcast;  heat  w’as  supplied  to  the  house  from  the  storage  v;all. 
(The  au:ciliary  furnace  is  in  the  basement  v;hich  stores  heat  that  is  not  in- 
cluded in  the  heat  balance  and  reduces  the  apparent  heat  load  on  this  day.) 

Figure  14  shovjs  the  hourly  temperature  variations  and  solar  insolation 
for  the  first  week  in  March,  There  v;ere  only  two  days  in  the  middle  of  the 
v/eek  that  v;ere  clear.  Figure  15  depicts  the  corresponding  thermal  behavior 
during  this  period.  Heat  stored  during  tvjo  clear  days  at  mid-week  is  suf- 
ficient to  carry  the  heat  load  for  the  next  three  partially  overcast  days 
because  the  ambient  temperatures  are  moderate. 

Note  that  the  largest  temperature  svjings  occur  on  very  cold,  clear  days. 
Since  heat  is  transferred  from  the  storage  wall  to  the  house  by  natural  con- 
vection, this  limit  on  the  rate  of  heat  transfer  results  in  a lov;er  minimum 
house  temperature.  The  higher  maximum  is  due  to  direct  sola^r  gain, 

10.0  PlLRFOrjvIANCL  SUiai/J^Y  DURING  MONITORING  PMRIOD 

Hourly  data  was  recorded  betvjeen  November  28  and  March  6,  A summary 
of  this  entire  period  is  shown  in  Table  3.  The  heating  degree  days  during 
the  test  period  ;jere  27  percent  above  average  creating  a challenging  test 
environment.  The  overall  average  house  temperature  vjas  19.0°C  ( 66°P)  vjith 
an  average  differential  of  3.7°C  (6.6°F)  betvieen  the  upper  and  lov/er  floors. 
The  average  daily  sizing  in  house  temperature  vzas  6.4°C  (11,5°F). 

During  the  test  period,  electrical  power  dissipated  in  the  house  was 
7.7  GJ  (8  X 10  BTU) . Thirteen  fires  v;ere  built  in  the  auxiliary  wood 
furnace  which  put  a total  of  10.3  GJ  (IO.9  x 10^  BTU)  into  the  house  (about 
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0,8  cord  of  pine  was  consvuned) , The  total  heat  requirement  of  the  house 
vias  52  GJ  ( 55  X 10^  BTU)  based  on  the  experimental  heat  loss  factor  of  O.83 
(486  BTUh”^ ) . Based  on  these  measurements,  the  sun  provided 
34  GJ  ( 36  X 10^  BTU)  or  65  percent  of  the  total  heating  energy  during  De- 
cember, January  and  February. 

The  average  efficiency  of  the  furnace  v;as  determined  by  weighing  the 

total  v/ood  put  in  during  each  firing  and  then  measuring  inlet  and  outlet 

duct  temperatures  and  outlet  flow  rate.  The  heating  value  of  the  wood  was 

—1  —1 

taken  as  I8.4  MJkg  (8,823  BTUlb  ).  Using  these  data,  the  average  overall 
efficiency  of  the  v:ood  furnace  vjas  58»7  percent  during  eleven  test  firings. 

11.0  PROJECTED  ANNUAL  PERFOPJIAhCK 

During  the  entire  heating  season,  a total  of  sixteen  fires  v/ere  built 
in  the  auxiliary  viood  furnace,  the  first  on  November  11  and  the  last  on 
January  27.  A total  of  11 40  kg  (25O8  lb,  about  one  cord)  of  wood  was  burned 
producing  about  12,7  GJ  (13.4  x 10  BTU)  of  heat.  The  annual  electric 
pov/er,  dissipated  as  heat  into  the  house,  contributed  about  15  GJ,  or  I6  x 
10  BTU  (700  ki-'h  mo  ).  Using  an  annual  value  of  5OOO  C days,  the  heat  load 
of  the  house  would  be  39,6  GJ  (105  x 10^  BTU)  and  the  annual  solar  fraction, 
72  percent.  This  summary  data  is  also  given  in  Table  3. 

12.0  SUI'fl/IARY 

The  test  results  shox^r  that  this  system  provides  a high  percentage  of 
the  house  heat  load  by  solar  conversion,  even  during  extremely  cold  weather. 
It  seems  unlikely  that  the  overall  solar  fraction  could  be  significantly  im- 
proved by  economically  admissible  changes  in  either  aperture  area  or  stor- 
age capacity.  The  conversion  efficiency  of  the  XT?ood  furnace  and  the  heat 
loss  factor  of  the  house  itself  could  not  be  greatly  improved. 

Reductions  of  over  10  percent  in  the  overall  heat  loss  vjould  result 
if  the  movable  insulation  could  be  constructed  so  as  to  have  a heat  transfer 
coefficient  equivalent  to  a normal  wall.  The  effective  insulating  value  of 
the  curtain  is  about  R— 5.  The  curtain  is  currently  opened  and  closed  by  a 
manual  sv;itch.  The  curtain  would  need  an  automatic  control  to  make  this 
system  acceptable  for  general  use.  The  furnace  is  not  automatic  and  this 
is  another  inconvenience.  A reduction  in  the  24-hour  temperature  variations 
would  increase  the  comfort  of  the  house.  This  is  a difficult  and  fundamental 
problem  of  passive  design. 
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Condensation  betv;een  the  curtain  and  the  window  is  a problem  and  mao'" 
eventually  cause  the  v;ood  sills  to  deteriorate.  The  insulating  curtain 
shovjs  signs  of  v;ear  and  v;ill  need  to  be  periodically  repaired  and  main- 


tained 
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TABLE  1 

CALCULATED  HEAT  LOAD,  PASSIVE  HOUSE 


R 

U 

Area 

U X A 

(Btu/hr  ft^ 

°P)  (sq.  ft.) 

Pioof 

40 

.025 

1112 

28 

Walls 

20 

.05 

1660 

83 

Floor 

40 

.025 

1144 

29 

VJ  indov;s 

1-4 

.68 

217 

148 

Skylights 

2.7 

.36 

30 

10.8 

Aperture: 

(Open) 

(1.4) 

(.68) 

(426) 

(290) 

Closed 

7.4 

.13 

426 

55 

Infiltration  (|-air 

change/hour  ~ 

• 20,592  ft^) 

20,592 

ft^  X -jVhr  X . 

018  = 

155 

TOTAL  HEAT 

LOAD: 

(With  Curtain  Open) 

744  Btu/hr°P 

(iJith  Curtain  Closed)  5^9 

Heighted  Average  571 


1.4  Mj/hr°C 
.96 


1.08 
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("b)  Interior  view  shov/ing  insulating  curtain 
and  storage  v/all 


Figure  1:  Fowlkes  Passive  House  Photographs 
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Figure  2:  Photographs  of  the  Passive  Storage  Wall  (c)  Solar  wall  as  seen  from  lower  floor  dressing  room. 
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(a)  Auxiliary  wood  furnace 


(b)  Hot  and 
cold  air 
ducting  to 
solar  wall 
(before 
tubes) 


Figure  3:  Auxiliary  Heating  System  Photographs 


Figure  L\  Floor  Plar  for  ground  Floo 
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Figure  5:  Floor  Plan  for  Second  Floor 
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Figure  6:  Shading  Diagram  Passive  House 
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Figure  7s  Section  Viev;  and  Transducer  Layout 
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Figure  8:  Solax  Wall  Transducer  Locations 
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PUT^ACE  OUTLET 


51  cm 


Average  Velocity 


5.63m/ 


sec. 


Volume  Flow  Rate 


.58lfTi^/sec  or  Volume  Plow  Rate  = 34»8m^/i 


min 


Figure  9t  Furnace  Outlet  Duct  Velobity  Distribution 


cm 


Curtain  down 
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Figure  10:  Hourly  Heat  Balance  Ec^uation 


lOOOREi'l  SUBR  f**  COMPUTE  HOURLT  SUMS 


3005  SUMMING" 

3070  N0=N0+1\R£M  HHRS 

3015  S(1)=D(T)t3.6\  REM  SOLAR  INSOLATION  INSIDE  MJ  PER  SM 
3020  S(2)=D(4)+D(1 )*127.0\REM  SOLAR-TUBE  ENERGY 
3030  S(3)=D(1)+36.0\R£M  SOLAR-HOUSE  ENERGY 
3040  S7  --(D(2)+l!(3))t3.6\REM  ELECT  UTIL  ENERGY 
3050  S<4)=S7  +0.6\R£M  ELECT  POUER  HEAT 

3060  S(5)=(D(7)-D(6) )*D(5)*2.41\REM  UOOD  FURNACE  HEAT, NEU  CALIBRATION 
3070  S(10)^(D(8)+2*Ii(10)  )/3\R£M  WEIGHTED  AVERAGE  TUBE  TEMP, GOOD  PROBES 
3075  S(11 )=D(19) 

3080  S(12)=(D(12)+D(13)+D(14)+D(15))/4.0\REM  HOUSE  TEMP  (AVG) 

3085  S(13)=(D(12)+D(13)-D(14)-D(15))/2  \REM  HOUSE  UP-DOUN  TEMPERATURE 
3090  S(6)=(S(10)-S8  )*52.0\REM  TUBE-STORED  HEAT, NEW  CONSTANT 

3100  S(7)  = (S(12)-S9  )+-8.0\REM  HOUSE-STORED  HEAT 

3110  S(8)=Sn2)-D(19)\R£M  DEGREE-HOURS 

3120  S(9.)  = (S(2)+S(3)+S(4)+S(5)-S(6)-S(7)  )/S(8)\REM  HEAT  LOSS  FACTOR 
3125  S(1 4)=S(10)-S(12)\REM  TUBE  MINUS  HOUSE  TEMP 
3130  *#1 ,23I,R,X4I,H0,Z3I,D0,X5I,H0,TAB(15), 

3140  FOR  1=1  TO  14\T(I)=T(I)+S(I)\REM  SUM  FOR  DAY,  PRINT  HR  DATA 
3150  '#1  ,%7F1 ,S(I),\NEXT\'»1 

3160  S8=S(10)\S9=S(12)\RETURN  \REM  SAVE  FOR  TEMP  DIFFERENCE  - TU9EXH0USE 
3200REM  CALC  AVERAGES  S FINAL  QUANTITIES 

3210  FOR  1=9  iO  14\  REM  AVERAGE  TEMPERATURES  (NOTE  HEATS  NOT  INCL.,) 

3220  T(I)=T(I)/NO\R£M  DAILY  AVGS  = SUM  / NR  HRS 
3230  NEXT\RETURN 
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3510  !#1\!#1  ,"DAY",%4I,M1  .3:31,01 ,3:5I,NO,TAB(15). 
3520  FOR  1=1  TO  14 

3530  !81,Z7F1.T(I),\NEXT\!tt1,"  SUM-AVG  "\'«1 
3540  RETURN  . 

3550END 


Figure  11:  Computer  Program  of  Energy  Ecjuations 
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Figure  12;  Hourly  Temperatures  During  A Heek  In  January 
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Figure  13j  Daily  Energy  Summary  Prom  Januaiy  4 Through  January  10 
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Figure  14:  Hourly  Temperatures  During  A Feek  In  March 
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Figure  15:  Daily  Energy  Suiranary  From  February  28  Through  March  6 
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Bozeman,  MT  59715 
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NOTICE 

This  report  was  prepared  as  an  account  of  work  sponsored  by  the  Energy 
Division  of  the  Montana  Department  of  Natural  Resources  and  Conservation 
through  the  Alternative  Renewable  Energy  Sources  Program.  Neither  the  State 
of  Montana,  nor  the  Department,  nor  any  of  their  employees,  nor  any  of  their 
contractors,  subcontractors,  or  their  employees,  make  any  warranty,  ex- 
pressed or  implied,  or  assxnnes  any  legal  liability  or  responsibility  for 
the  accuracy,  completeness  or  usefulness  of  any  information,  apparatus,  pro- 
duct or  process  disclosed,  or  represents  that  its  use  would  not  infringe 
privately-owned  rights. 
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This  project  is  an  active-liquid  solar  system  retrofit  Belonging  to  Dr. 
Richard  Sheridan,  Missoula,  Montana.  The  free-stamding  collectors  are  lo- 
cated in  a meadow  about  200  feet  from  the  house.  The  collectors  deliver 
heat  to  water-filled  heat  storage  tanks.  A second  liquid  loop  carries  stored 
heat  to  an  exchanger  in  the  air  duct  which  in  turn  delivers  heat  to  the  house. 

The  solar  system  does  not  deliver  solar  heat  in  mid-winter  due,  in  part, 
to  shading  of  the  collector  array.  During  the  spring  the  collectors  de- 
livered 20^  of  the  available  solar  insolation  to  the  heat  storage  but  only 
3.6^  to  the  house.  The  heat  lost  from  the  storage  system  did  not  go  into  the 
house.  There  appeared  to  be  design  and  control  problems  in  the  storage-to- 
house  heat  exchange  loop.  The  solar  system  supplied  of  the  heat  load  of 
the  house  during  the  monitoring  period.  It  is  estimated  that  system  changes 
could  increase  the  solar  fraction  to  25^  on  an  annual  basis. 

The  solar  domestic  hot  water  system  was  also  monitored.  It  was  a much 
better  performer  that  the  solar  heating  system.  The  DHW  solar  system  de- 
livered an  average  of  28^  of  the  solso*  insolation  to  the  preheat  tank.  This 
amounts  to  an  average  of  15  mega  joules  ( 16,000  Btu)  each  day. 


NOTE  ON  UNITS 

The  data  is  presented  in  metric  imits.  These  units  are  in  accord  with 
the  recommendations  of  the  International  Solar  Energy  Society  as  well  as 
most  current  technical  journals.  The  overall  summaries  and  numbers  in  the 
test  have  been  given  in  dual  units.  The  heat  load  calculations  are  pre- 
sented in  conventional  units. 

Most  readers  are  familiar  with  the  Celsius  or  Centigrade  temperature 
scale  and  aj'e  also  familiar  with  energy  \inits  of  kilowatt  hours.  The  only 
other  metric  unit  which  occurs  frequently  in  this  report  is  the  Joule  ( j) , 
the  metric  unit  of  energy  which  replaces  the  Btu  (a  Btu  is  equal  to  1,055 

Joules).  The  mega  joule  (MJ)  is  equal  to  lO^J  (100,000  Btu/hr  = 105  Mj/hr). 

2 

Other  abbreviations  used  in  this  report  include:  m = meters;  m = 
square  meters;  1 = liters;  h = hours;  gph  = gallons  per  hour;  and  MBtu  = 
millions  of  Btu*s. 
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1.0  INTRODUCTION 

This  project  is  in  the  residence  of  Dr,  Richard  Sheridan  and  is  lo- 
cated on  the  north-facing  slope  of  Pattee  Canyon  near  Missoula.  Dr.  Sheri- 
dan, in  addition  to  being  a Professor  at  the  University  of  Montana,  is  an 
active  environmentalist  and  solar  advocate.  He  is  also  an  energetic  and 
skilled  craftsman  and  has  designed  and  built  his  own  house  as  well  as  the 
solar  heating  system  monitored  in  this  report. 

The  house  was  sited  and  built  over  10  years  ago.  Solar  orientation  and 
energy  conservation  were  not  among  the  primary  design  considerations  at  that 
time.  When  Dr.  Sheridan  decided  to  retrofit  solar  collectors  to  his  house, 
he  located  the  collectors  in  a meadow  about  200  feet  from  the  house  in  order 
to  gain  exposure  to  the  sun, 

2.0  ACKNOWLEDGMENT 

The  cooperation  and  advice  of  Dr.  Sheridan  was  crucial  to  the  success 
of  this  monitoring  project.  He  kept  records  of  the  readings  of  the  manual 
flow  meters  which  were  necessary  to  characterize  the  average  flow  rates  of 
the  pumps.  He  contributed  several  suggestions  concerning  the  behavior  of 
the  system  and  helped  with  the  practical  problems  of  installing  the  probes 
and  routing  the  extension  wires.  The  data  tapes  were  changed  by  Dr,  Sheri- 
dan and  he  alerted  us  when  there  were  problems.  We  gratefully  acknowledge 
his  contributions  and  enthusiastic  support  of  this  monitoring  project, 

3.0  DESCRIPTION  OF  SOLAR  SYSTEM 

A schematic  of  the  solar  system  is  shown  in  Figure  1 and  a house  plan 

is  shown  in  Figure  2,  The  main  collector  axray  consists  of  24  licjuid- 
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cooled  panels  totalling  37*9  (408  ft  ) in  area  and  tilted  at  45 

facing  south.  Figure  3.  The  collectors  are  double  glazed  and  have  an  insu- 
lated wood  frame.  The  absorber  plate  is  painted  flat  black  and  is  attached 
to  copper  tubes  about  10  cm  (4  in)  on  center.  Heat  transfer  putty  was  used 
between  the  tubes  and  the  absorber  plate.  This  material  has  been  unsatisfac- 
tory as  it  runs  out  of  the  joint  at  higher  temperatures  as  shown  in  Figure 
4(a).  The  absorber  plates  are  manifolded  and  the  fluid  carried  to  the  stor- 
age tanks  in  the  house  through  pipes  enclosed  in  foam  insulation.  Figure  3(d), 

The  collector  fluid  is  a 505^  ethylene  glycol-water  mixture  having  a den- 
sity of  1.07  kg/l  and  a specific  heat  of  3.18  kj/kg°C.  The  collector  pump 
is  powered  by  a 3/4-hp  electric  motor  and  produces  a flow  of  2089  Ipb  (552 


2 


gph) . The  collector  fluid  passes  through  finned  tubes  immersed  in  three 
fiberglass  heat  storage  tanks  each  filled  with  5»677  1 (500  gal)  of  water* 

One  tank  is  divided  so  that  storage  mass  can  be  adjusted  to  four  levels  by- 
setting  appropriate  values,  Figure  1.  The  heat  storage  tanks  are  located 
in  a basement  garage  and  are  individually  insulated  with  R-11  fiberglass. 
Figure  5 shows  the  tanks  and  the  circulating  pumps. 

Heat  removal  from  storage  is  shown  in  the  schematic  diagram  in  Figure 
1.  Fluid  is  circulated  by  a l/l2-hp  pump  which  delivers  a flow  of  5^7  Iph 
( 1 55  gph) . The  fluid  passes  through  finned  tubes  in  the  heat  storage  tanks 
and  then  through  finned  tubes  located  in  an  air  duct  supplying  the  house, 
Figure  5(a)*  An  auxiliary  hot-air,  oil  furnace  is  also  connected  to  this 
duct.  Additional  auxiliary  heat  is  supplied  by  a wood-burning  stove.  Fig- 
ure 1 . 

The  domestic  hot  water  solar  pre-heat  system  is  also  shown  in  Figure  1. 

The  collectors.  Figure  4(b),  are  similar  in  construction  to  the  main  collec- 
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tor  array  and  have  an  area  of  4*7  m (51  ft  ).  A 50?^  ethylene  glycol-water 
mixture  is  circulated  by  a pump  delivering  an  average  flow  of  250  Iph  (66 
gph) . This  flow  is  adjusted  manually  by  the  owner  to  control  the  collector 
outlet  temperature.  The  collector  fluid  passes  through  an  exchanger  in  a 
pre-heat  tank  having  a capacity  of  405  1 (105  gal).  VJell  water  enters  the 
pre-heat  tank  and  then  goes  through  a conventional  electric  water  heater. 

3. 1 Monitoring  Instrumentation 

The  transducer  arrangement  is  shown  in  the  system  schematic  diagram  in 
Figure  1 . This  solar  system  has  many  components  and  heat  transfer  loops 
making  characterization  and  instrumentation  a formidable  problem.  The  owner 
was  interested  in  the  performance  of  the  solar  domestic  hot  water  system 
and  believed  that  it  would  prove  to  be  more  efficient  than  his  solar  space 
heating  system.  Because  of  these  considerations,  the  available  instrumen- 
tation resources  were  focused  on  the  solar  components  of  the  system;  the  oil 
furnace  and  wood  stove  heating  systems  received  secondary  attention. 

The  solar  collector  loop  was  instrumented  with  a continuous  reading  flow 
meter,  a status  relay  on  the  pump  motor  and  three  temperature  transducers. 

Two  of  these  transducers  measured  inlet  and  outlet  collector  temperatures  at 
the  storage  tank  to  determine  solar  heat  delivered  to  the  storage.  An  extra 
probe  measured  the  collector  outlet  temperature  at  the  point  where  the  total 
collector  flow  left  the  array  and  started  toward  the  house,  location  I)  in 
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Figure  3.  This  allowed  an  estimate  of  the  heat  loss  due  to  the  remote  lo- 
cation of  the  collector  array.  Solar  insolation  was  measured  at  the  center 
of  the  array  and  in  the  plane  of  the  collectors. 

The  average  temperature  of  the  thermal  storage  was  sampled  with  three 
probes  attached  to  the  surfaces  of  the  three  heat  storage  tanks.  The  probes 
were  covered  with  heat  transfer  putty  and  insulated.  Net  heat  changes  in 
the  storage  system  are  reflected  by  the  change  of  these  temperatures  with 
time. 

The  heat  extraction  loop  was  instrumented  with  a manual -reading  flow 
meter,  a status  relay  on  the  pump  motor  and  temperature  probes  on  the  inlet 
and  outlet  pipes.  The  DHW  collector  loop  was  similarly  instrumented.  A 
probe  on  the  DHW  storage  tank  reflected  energy  changes  in  the  tank.  Elec- 
tric power  is  monitored  continuously  as  well  as  house  temperature  and  am- 
bient temperature.  A probe  in  the  furnace  outlet  duct  allowed  for  an  ap- 
proximate calculation  of  auxiliary  furnace  heat  output. 

4.0  DATA  BASE  AND  ACCURACY 

The  relative  accuracy  of  the  temperature  measurements  was  0.2®C  and  the 
accuracy  of  temperature  changes  was  better  than  O.I^C.  The  flow  meters  are 
accurate  to  within  3^.  The  electric  power  measurements  are  calibrated 
against  the  watt  meter  and  are  correct  within  3/^.  The  solar  insolation  is 
accurate  to  within  on  the  basis  of  average  daily  total  insolation.  The 
accuracy  of  the  calculated  heat  flows  changes  depending  on  the  size  of  the 
temperature  difference.  On  the  average  the  error  is  believed  to  be  no  more 
than  15/0. 

The  data  acquisition  system  was  installed  March  19»  1979  and.  removed 
June  19,  1979»  A total  of  1919  hours  of  data  were  recorded  for  I8  trans- 
ducer channels  of  34,542  hourly  average  data  values.  Each  hourly  average 
consisted  of  about  66O  individual  readings;  hence,  the  data  base  contains 
about  23,000,000  channel  readings.  The  total  missing  data  was  12^  and  was 
primarily  due  to  a data  tape  that  was  not  replaced  when  it  ran  out. 

5.0  MONITORING  DATA  AND  RESULTS 

5. 1 Site  Survey 

A shading  profile  taken  in  the  center  of  the  collector  array  is  shown 
in  Figure  5«  The  arcs  show  the  path  of  the  sun  during  winter  between  the 
fall  and  spring  equinox  for  various  times  of  day.  The  composite  photograph 
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shows  the  view  of  the  sky  looking  south,  seen  hy  the  collector  axray.  The 
shading  profile  indicates  that  there  is  a high  percentage  of  shading  during 
the  months  of  November,  December  and  January.  The  implications  of  this 
shading  diagram  are  substantiated  by  Dr,  Sheridan’s  reports  that  the  solar 
system  puts  out  little  or  no  heat  during  this  season. 

5.2  House  Heat  Load 

The  calculated  heat  load  factor  of  the  house  is  shown  in  Table  1. 
Dominating  features  affecting  the  heat  load  include  a large  cathedral  win- 
dow on  the  north  end  of  the  building  and  a large  roof  area  inherent  in  the 
A-frame  design,  which  is  insulated  to  R-10.  Dr.  Sheridan  reports  that  he 
reduces  the  heat  requirements  of  the  house  by  (a)  closing  off  bedrooms, 

(b)  reducing  the  thermostat  when  the  house  is  unoccupied  and  at  night,  and 

(c)  keeping  the  house  temperatures  around  16°C  (60°F).  The  data  demonstrates 
some  of  these  patterns, 

5.3  Flow 

The  manual  flow  meter  readings  taken  daily  by  Dr,  Sheridan  axe  repro- 
duced in  Table  2,  The  net  daily  flows  were  divided  ly  the  daily  on-time  of 
the  respective  pumps  to  yield  a daily  average  flow  rate  for  each  system,  A 
large  number  of  these  readings  were  averaged  to  produce  flow  rate  values 
which  were  subsequently  used  to  perform  the  hourly  heat  balances.  Checks  on 
the  heat  balance  of  the  storage  system  showed  these  values  to  be  reliable. 

The  flow  meters  have  an  accuracy  of  better  than  3^  and  their  calibration  was 
checked  after  removal  from  the  test  site.  The  calibration  facility  at  the 
Bozeman  City  Water  Department  was  used  for  these  determinations, 

5«4  Collector  Efficiency 

Efficiency  curves  for  the  Sheridan  collectors  are  shown  in  Figure  6 
along  with  a HUD  reference  curve  for  this  general  type  of  collector.  The 
”output’'  of  the  test  collector  was  defined  as  heat  delivered  to  the  storage 
system.  The  collector  is  thus  penalized  for  heat  losses  in  the  extensive 
piping  and  manifold  system  between  the  collectors  and  the  storage  tanks. 

This  loss  appears  to  be  around  10^  of  the  output.  The  temperature  rise  of 
the  liquid  flows  through  the  collector  was  typically  5°  to  10°C  (9°  to  18°F). 
The  temperature  drop  between  the  collector  array  outlet  and  the  storage  tank 
inlet  was  around  0.5°C  (l°F).  This  indicates  about  a 5^  loss  in  that  section 
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of  piping.  The  losses  in  the  manifold  piping  within  the  collector  field 
appeax  to  he  another  5^. 

5.  5 Heat  Exchanger  Effectiveness 

Table  3 shows  a modified  version  of  the  final}  hourly  sximmary  program 
for  two  typical  days.  The  columns  of  interest  in  this  table  are  noted  (I)} 

(2)  and  (3)  and  give  information  on  the  effectiveness  of  the  three  primary 
heat  exchangers.  The  value  of  the  Z\t  (temperature  differences)  leaving 
each  of  the  heat  exchangers.  If  the  heat  exchangers  were  ^OOfo  effective, 
the  fluid  leaving  the  exchanger  would  be  equal  to  the  heat  exchanger  tem- 
perature and  these  Z\t  values  would  be  *’0". 

The  exchangers  (1)  and  (2),  ’’into”  and  "out  of"  the  heat  storage  tank 
are  finned  tube  liquid  exchangers.  For  the  collector  exchanger  (I),  the 
temperature  change  through  the  exchanger  (not  shown)  is  of  the  same  order 
as  the  Z\T  leaving  the  exchanger.  For  exchanger  (2),  the  temperature 
change  through  the  exchanger  (not  shown)  is  around  ,25  to  ,5  of  the  /\ T 
leaving  the  exchanger.  The  other  end  of  exchanger  (2)  is  a finned  tube  in 
the  house  air  duct  ( 3) . The  Z\t  between  the  exit  fluid  and  the  air  is 
around  20°C  ( 36°F) . 

These  data  indicate  that  the  heat  exchange  loop  from  the  storage  to 
the  house  has  a low  effectiveness  at  both  ends  and  is  not  doing  an  efficient 
job  of  moving  the  heat  from  the  storage  to  the  house.  The  collector-to- 
storage  heat  exchanger  does  a much  more  efficient  job. 

5.6  Hourly  Performance  Sivnmary 

The  fundamental  unit  of  the  performance  data  is  the  hourly  summaiy.  This 
summary  includes  all  the  individual  heat  flows  and  allows  determination  of 
their  contributions  to  the  overall  system.  Figure  8 lists  the  portion  of  the 
computer  program  used  to  process  the  hourly  data  to  compute  heat  quantities. 

Figure  9 explains  the  quantities  tabulated  in  the  overall  performance 
summary.  In  general,  this  summary  traces  the  hourly  average  flow  of  energy 
from  solar  radiation,  to  collectors,  to  storage  and  to  the  house  (or  the  DHW 
tank).  Auxiliary  heat  inputs  due  to  electrical  dissipation  and  the  furnace 
are  tabulated.  The  hourly  house  heat  loss  is  based  on  the  heat  load  factor. 
Table  1 , and  the  measured  difference  in  temperature  between  the  inside  and 
outside  of  the  house.  Temperatures  at  important  points  in  the  system  are 
also  tabulated. 
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The  hourly  performance  for  35  days  is  given  in  the  tables  in  Appendix 
I.  As  an  example  of  the  use  of  this  information,  look  at  Table  4 for  the 
day  of  April  24.  The  sun  came  up  between  7 and  8 AM  and  set  after  8 PM ; 
see  ’’SOLAR  INSO"  column.  The  collector  pump  turned  on  shortly  before  11  AI^l 
and  delivered  heat  to  storage  until  4 PM  (18  hours);  see  ’’COLL  OUTPT”  col- 
umn. The  storage  output  column,  ’’STORE  OUTPT”,  shows  a loss  of  2.1  Mj/hr 
during  the  early  morning  hours.  The  storage  begins  delivering  heat  to  the 
house  between  10  and  1 1 AM  and  begins  receiving  significant  heat  from  the 
collectors  between  11  AM  and  noon. 

It  should  be  noted  that  the  hourly  "COLL  OUTPUT”  and  "STORE  HOUSE” 
quantities  are  much  more  accurate  than  the  "STORE  OUTPT”  quantities.  The 
storage  output  is  calculated  from  the  temperatures  measured  at  only  three 
points  in  the  storage  system.  The  temperature  differences  are  also  small 
and  one-half  hour  out  of  phase.  The  combination  of  these  features  causes 
the  hourly  values  to  lag  and  to  fluctuate.  The  average  values  taken  over 
one  or  two  days,  however,  balance  as  they  should:  i.e.,  Input  = Output  + 
Loss.  The  use  of  the  "STORE  OUTPT”  column  is  to  give  a rough  check  on  the 
"COLL  OUTPT”  AND  "STORE  HOUSE”  quantities  and  to  determine  the  average  heat 
loss  factor  for  the  heat  storage  system. 

The  "DHWC  OUTPT”  is  the  heat  delivered  by  the  domestic  hot  water  col- 
lectors to  the  heat  storage  tank.  "DELTA  DHW”  is  again  calculated  on  the 
temperature  change  of  one  probe  on  the  DHW  tank  and  shows  the  lag  and 
hourly  fluctuation  previously  discussed. 

The  "DELTA  DHW”  column  shows  that  the  heat  loss  was  less  than  0.05  Mj/hr 
during  the  early  morning  hours.  Between  7 and  9 AM  some  heat  was  removed 
from  the  tank  probably  due  to  normal  household  uses.  Between  10  and  1 1 AM 
the  DHW  collectors  turned  on  and  began  delivering  heat  to  the  tank.  Small 
amounts  of  heat  were  withdrawn  during  the  evening.  (Also  note  the  increased 
hourly  heat  loss  due  to  the  higher  tank  temperature.)  Between  10  PM  and 
midnight  (23  and  0 hrs.)  it  appears  that  there  were  significant  withdrawals 
from  the  tank,  perhaps  due  to  baths  or  washing  dishes  or  clothes.  Note  that 
the  electric  power  between  10  and  11  PM  (23  hrs.)  increased  markedly,  proba- 
bly due  to  the  auxilieny  DHW  heater.  The  "PURNC  DUCT”  column  shows  that  the 
furnace  was  used  for  about  3 hours  in  the  morning  when  the  Sheridans  got  up. 

The  hourly  electric  power  use  (divide  MJ  by  3.6  to  get  kWh)  is  typical 
in  that  it  shows  the  morning  and  evening  peaks.  The  steady  daytime  use  of 
about  2.5  MJ  (.7  Id/)  is  mostly  due  to  the  solar  circulating  pumps. 
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The  bottom  line  of  the  summary  shows  the  daily  sum  of  all  energy  items 
and  the  average  of  the  temperature  items.  The  total  heat  inputs  to  the  house 
on  this  day  consisted  of  112  MJ  from  solar  storage,  ^2  MJ  electric  power  and 
168  MJ  from  the  auxiliary  furnace.  This  total  of  332  MJ  compares  tolerably 
well  with  the  calculated  heat  loss  of  416  MJ.  It  is  possible  that  the  wood 
stove  was  being  used  in  the  evening  as  the  house  temperature  peaks  at  21.2°C 
(70°F)  at  9 PM. 

There  is  a wealth  of  information  concerning  the  dynamic  behavior  of  the 
solar  system  contained  in  the  hourly  and  daily  summaries  in  Appendix  I.  In 
order  to  summarize  this  information  all  the  data  in  Appendix  I have  been 
condensed  into  daily  summaries  covering  the  entire  monitoring  period.  These 
data  are  given  in  Tables  5,  6 and  7.  The  headings  and  units  are  the  same  as 
discussed  previously. 

6.0  OVERALL  PERFORI^IANCE  SmmARY 

This  data  has  been  further  summarized  to  produce  the  results  in  Tables 
8 and  Table  8 siammarizes  the  data  for  the  solar  heating  system  for  78 
data  days  between  March  20  and  June  l8,  1979»  This  data  shows  that  the  col- 
lectors delivered  20^  of  the  available  radiation  to  the  storage  tanks.  Of 
this  stored  solar  heat,  l8^  was  delivered  to  the  house,  the  remainder  being 
lost  to  the  garage. 

The  second  part  of  Table  8 compares  the  overall  heat  load  of  the  house 
to  the  relative  contributions  of  the  solar  system  and  the  auxiliary  system. 
This  table  shows  that  the  solar  heating  system  contributed  "jfo  and  the  auxil- 
iary systems  ( furnace  and  wood)  contributed  89^.  The  remaining  heat  was 
electrical  dissipation  due  to  appliances  and  lights,  estimated  as  0.6  of  the 
total  electrical  consumption. 

Table  9 summarizes  the  performance  of  the  domestic  hot  water  solar  sys- 
tem. Solar  insolation  was  not  measured  at  the  DHW  collector  panel.  Because 
of  its  shaded  location,  it  was  assumed  to  be  exposed  to  0,75  of  the  solar 
radiation  measured  in  the  main  collector  field.  Based  on  this  estimate, 
the  collector  delivered  28^  of  the  incident  radiation  to  the  tank.  If  the 
average  heat  loss  of  the  tank  is  considered,  the  useful  heat  output  amounts 
to  2/!ffo  of  the  incident  radiation. 
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7.0  COMMENTS  AND  CONCLUSIONS 

This  solar  heating  system  has  three  characteristics  which  reduce  its 
performance: 

(1)  The  collectors  are  shaded  in  mid-winter. 

(2)  There  are  heat  losses  of  about  10^  in  the  piping  between  the  col- 
lector array  and  the  house. 

(3)  Only  ^&fo  of  the  heat  that  does  get  into  the  storage  tanks  is  de- 
livered to  the  house. 

Items  (1)  and  (2)  could  not  be  avoided  due  to  the  sitting  of  the  house. 
Item  (3)  appears  to  have  three  main  causes: 

(a)  Heat  lost  from  the  storage  tanks  does  not  enter  the  house  envelope 
and  is,  thus,  truly  lost. 

(b)  The  heat  exchangers  between  the  storage  tanks  and  the  house  have 
limited  ability  to  transfer  heat  to  the  house.  Data  in  Table  3 as 
well  as  the  hourly  performance  suggests  that  these  exchangers  should 
be  of  higher  efficiency. 

(c)  The  storage-to-house  exchangers  do  not  operate  on  a large  percentage 
of  the  days  even  though  the  storage  is  charged  to  temperatures  of 
40-50°C  0 04-122°P);  see  Tables  5f  6 and  7. 

If  all  the  heat  collected  during  the  monitoring  period  could  have  been 
delivered  to  the  house,  the  solar  fraction  would  have  been  27^  instead  of  75^. 

The  domestic  hot  water  system  has  a better  overall  performance  in  terms 
of  useful  heat  delivered.  Since  the  storage  tank  is  located  inside  the  house, 
losses  from  this  tank  contribute  useful  heat  to  the  house. 

8.0  PREDICTIONS  OF  ANNUAL  PERFORMANCE 

In  this  section  the  performance  during  the  monitoring  period  is  pro- 
jected to  create  an  estimate  of  the  annual  performance.  The  starting  point 
of  this  analysis  is  the  environmental  data  concerning  the  solar  radiation 
available  and  the  ambient  temperatures  or  degree  days. 

The  solar  radiation  measured  at  the  Hellgate  High  School  in  Missoula  is 
compared  with  the  radiation  measured  at  the  Sheridan  site  in  Table  10.  The 
Missoula  data  is  taken  from  measurements  made  as  part  of  the  SBIM  program 
sponsored  by  the  Department  of  Natural  Resources  and  Conservation.  The  mea- 
surements were  taken  on  a 60°  tilted  surface.  This  data  is  mathematically 
adjusted  to  correspond  to  a 45°  tilted  surface,  the  plane  of  the  monitoring 
devices  on  the  Sheridan  collectors.  The  predicted  data  and  the  data  measured 
on-site  show  satisfactory  agreement. 
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Figure  10  is  a graph  of  the  performance  of  the  Sheridan  solar  collectors 
during  the  monitoring  period.  The  total  insolation  striking  the  collector 
field  is  plotted  against  the  heat  delivered  by  the  collectors  to  the  storage 
system,  A single  curve  has  been  drawn  through  these  points.  This  curve  is 
assumed  to  characterize  the  daily  output  of  the  collectors  given  a solar  in- 
put. The  scatter  of  the  data  points  shov/s  that  the  straight  line  is  only  an 
approximation  to  the  data. 

Table  11  compiles  the  projected  annual  performance.  Long-term  average 
degree  day/month  data  for  Missoula  is  listed  in  the  first  column.  This  data 
is  based  on  18,3°C  (65°F).  The  second  column  is  the  degree  day  data  adjusted 
for  a house  temperature  of  l6*^C  (60°P)  to  better  reflect  the  real  average 
temperature  of  the  Sheridan  house. 

The  next  column,  INSOL  @ 60°,  is  the  average  daily  solar  radiation  re- 
corded in  Missoula  on  a 60°  tilt  during  the  time  period  from  July  1978  to 
June  1979»  The  factor  in  the  next  column  adjusts  this  data  to  the  45°  tilt 
of  the  Sheridan  collectors.  This  value  multiplied  by  the  anea  of  Sheridan's 
collector  array  and  converted  to  mega  joules  is  listed  in  the  next  column, 
TOTAL  INSO.  This  is  the  total  radiation  intercepted  by  the  collectors  as- 
suming no  shading. 

This  total  radiation  is  then  converted  to  the  collected  radiation  using 
Figure  10  and  is  listed  as  SOLAR  COLLECT  in  the  next  coliimn.  The  next  col- 
umn is  the  fraction  of  the  collected  energy  that  is  actually  delivered  to  the 
house,  based  on  the  average  performance  during  the  monitoring  period.  The 
heat  load  is  then  calculated  for  the  average  day  based  on  the  heat  load  fac- 
tor in  Table  1 and  the  adjusted  degree  days.  If  the  solan  collected  or  de- 
livered is  greater  than  the  heat  load,  the  excess  heat  is  not  usable.  In 
these  cases  the  heat  load  is  shown  in  parentheses  in  the  solar  col\xmn  and 
represents  the  maximiun  usable  heat  delivered.  The  solar  fraction  on  these 
days  is  1,00  or  100^^  and  can  be  no  greater. 

This  analysis  predicts  that  the  annual  solar  heat  delivered  by  the  sys- 
tem v;ill  be  5»760  MJ  (6,2  X 10^  Btu)  or  6^  of  the  annual  heat  load.  The 
analysis  shows  that  insignificant  amounts  of  heat  are  delivered  during  No- 
vember, December  and  January,  These  estimates  axe  in  general  agreement  with 
Dr.  Sheridan's  reports  on  the  actual  output  of  the  system  during  these  months. 

If  the  solar  heat  that  is  currently  being  delivered  to  the  storage  sys- 
tem, SOLAR  DELVD,  could  be  effectively  used  in  the  house,  the  solar  fraction 
would  be  much  larger.  This  assumption  is  represented  in  Table  11  for 
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comparison.  These  results  predict  that  2^o  of  the  heating  load  would  he 
met  by  the  soIclt  system  in  this  hypothetical  situation.  This  could  be  prac- 
tically accomplished  by  locating  the  heat  storage  inside  the  house  instead 
of  the  garage  and  increasing  the  effectiveness  of  delivering  stored  solar 
heat  to  the  house. 

9.0  GRAPHICAL  PRESENTATION  OF  RESULTS 

A graphical  presentation  of  the  hourly  data  for  two-five  day  periods  is 
presented  in  Figures  11  and  12.  The  figures  contain  data  for  solar  insola- 
tion and  four  temperatures.  The  general  pattern  shows  the  daily  increase  in 
storage  tank  temperatures  due  to  the  addition  of  solar  heat.  The  decrease 
in  storage  temperature  reflects  heat  loss  or  an  intentional  withdrav7al  of 
heat.  The  slopes  of  the  storage  temperature  at  night  show  that  the  DHl^  tank 
loss  rate  is  less  than  the  main  storage  loss  rate.  Note  the  small  tempera- 
ture drops  during  the  mornings  of  4/28  and  4/29  in  Figure  11.  These  corres- 
pond to  heat  withdrawals  by  the  house  pump.  The  columns  STORE  HOUSE  in 
Tables  4 and  5 in  Appendix  I show  the  amount  of  heat  taken  from  storage  (as 
measured  in  the  storage-to-house  heat  exchanger) . 

Note  that  the  DHW  tank  temperature  shows  large  energy  withdrawals  most 
mornings  and  evenings.  This  is  a typical  household  use  pattern.  The  graphs 
illustrate  that  the  solar  DHW  heat  is  being  "used"  while  most  of  the  solar 
heating  energy  is  being  dissipated  gradually  into  the  garage.  The  data  ta- 
bles (Appendix  l)  give  more  quantitative  results  for  the  energy  flows. 

Graphs  such  as  those  shown  in  Figures  11  and  12  are  useful  for  de-bugging 
and  checking  the  data.  They  are  also  fun  to  look  at. 
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TABLE  1 


HEAT  LOAD  OF  SHERIDAN  HOUSE 

R 

U AREA 

U X A 

(Btu/hr  ft^  °F)  (sq.  ft.) 

Roof 

10 

.10 

1757 

176 

Walls 

15.2 

.065 

1224 

79.5 

*Floor 

30 

.033 

1288 

43 

Windows 

1.5 

.65 

324 

210.6 

♦Infiltration 

i:  18,676 

ft^  X 1 X .018 

168 

Total  Heat  Load  677  Btu/hr  °F 

or 

1.28  Mj/hr  °C 


^Effective  R due  to  insulated  iDasement 
air  change/hr 


DATE  I.IETER(I)  METER(2)  ™w(1)  FL0W(2)  — TABLE  2 — DATE  METER(1)  METER(2)  PL0W(1)  FL0W(2) 

MAR  22  1 15.00  54.00  700.00  310.00  i MAY  6 1 1 12.00  912.00  120.00  10.00 
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TABLE  8 

SOLAR  HEATING  SIBIMARY  FOR  78  DATA  DAYS,  3/20  TO  6/l8 


MJ 

MBTU 

Total  Solar  Insolation 

44,663 

47.1 

100 

Solar  Collected  and  Stored 

9,077 

9.58 

20 

Solar  to  House 

1,628 

1.72 

3.6 

Lost  from  Storage 

7,449 

7.86 

House  Heat  Load 

24,252 

25.59 

100 

( 1 ) Furnace  Heat 

14,302 

15.09 

59 

(2)  Electric  Dissipation 

1,806 

1.91 

8 

Solar  Heat  Delivered 

1,628 

1.72 

7 

(3)  VJood 

6,516 

6.87 

26 

Average  House  Temperature 

• 

CO 

65.3°F 

(1)  Approximated  (-18^) 

(2)  ABsmning  ,6  of  total  electricity  dissipated  as  heat  inside  house 

(3)  Estimated  from  heat  balance  considerations 
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TABLE  9 

SOLAR  HOT  WATER  SUT^IARY  FOR  78  DATA  DAYS 


(1)  Total  Dffi'/  Insolation 

4218  Mj/m^ 

(1) 

Estimated  from  insola- 
tion measured  at  col- 

Total Collector  Output 

1178  MJ 

lector  array. 

(2)  DHW  Tank  Loss  (19^) 

177  MJ 

(2) 

Heat  "loss"  goes  into 
insulated  basement 

(3)  Daily  Average  Output 

15  MJ 

(3) 

15  Mj/day  is  enough  to 

Average  Collection 

2&fo 

heat  20.4  gallons  each 

Efficiency 

day  from  'to  135°P 

Average  System  Efficiency 

24/0 

TABLE  10 

PREDICTED  SOLAR  RADIATION  AT  SHERIDAN  SITE,  k^Jh/m^-day 


April 

May 

(2)  Jvine 


(1)  Missoula 
@ 60°  Data 

45°/60*^  Factor 
ASHRAE 

Predicted 

@ 45° 

Measured  @ 45* 
On  Site 

3.59 

1.07 

3.84 

3.65 

4.15 

1.11 

4.60 

4.79 

4.55 

1.13 

5.14 

4.69 

(1)  Montana  Solar  Data  Manual,  198O  edition 

(2)  Two  weeks  data 


9-? 


Figure  ?:  Sketch  of  Sheridan  residence  showing  window  areas 


?3 


(a)  House 


(B)  Collector  Array 

(C)  Distribution  Manifold 

(D)  Collector  Flow  - To  and 
from  house. 


Figure  3;  Photographs  of  main  collector  array 
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(a)  CloseMip  of  collector  panel 


(b)  Domestic  hot  water  (DHW)  collector 


Figure  4:  Collector  panels 
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(c)  Insulated  heat  storage  tanks 


Figure  5j  Heat  storage  system  and  controls,  located  in  basement 
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Figure  6;  Shading  diagram  showing  exposure  of  collector  field  to  the  sun  during  the  winter 
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$ *iOMI  01.4,43 


Figure  "J:  Collector  efficiency  curves 
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;?000R!-:f1  SUBR  + + COHPUTfc'  HOURLY  SUMS 
3005  SUMMING" 

3010  N0  = N0+1\REH  ItHRS 

3012  REM  ALL  THE  FO! LOUING  ARE  IN  MEGA  JOULES 

3015  S(1 ) =0(2) +.97+37. 9+3. 6 \R£M  SOLAR  INPUT 

3018  REM  TEMPERATURES  ZEROED  ON  STATUS,  GLYCOL  MIXTURE 

3020  S(2)=D(5)  + 2. 089+3. -l+(D(20)-D(16)  )/(D(5)  + . 0001  ) \ REM  COLL  TO  STORAGE 

3030  E(3)=5. 677  + 4. 19  + (D(12)-T1 ) \REM  MAIN  STORAGE  HEAT  CHANGE 
3040  S(4)=Ii(6)  + 4.  19  + 0. 587+(Ii(15)-D(19))/(Ii(6)  + . 0001  ) \REM  STORAGE  TO  HOUSE 
3050  S(5)=D(7)+0. 250+3. 4+(D{l7)-P(13) )/(D(7)+. 0001 ) \REM  DHU  COLLECTORS  OUT 
3060  S(6)=0. 405  + 4. 19+(D(18)-T2)  \REM  HOT  UATER  TANK  HEAT  CHANGE 

3070  S(7)=(D{3)+D(4))+3.6+.114  \REM  ELECTRICAL  POWER,  TOTAL  CORRECTED 

3071  S(8)-=0\IF(Ii(11  )-D(10))>0THEN  3072  ELSE  3073 


3072  S(8)=2+(Ii(11  )-D(10)  ) + .7\REM  ESTIMATED  AIR  DUCT  OUTPUT  .CORRECTED 

3073  09-D(12)-([i(9)+0.5  + D(10))/r.5 

3074  S(13)=-S(3)/'v'9  \ REM  STORAGE  LOSS  FAC  FOR 

3076  S(9)  = 1 .28+(D(10)-Ei(9) ) \REM  CALCULATED  HOUSE  HEAT  LOAD 

3080  T1=D(12)\T2=D(18)  \ REM  SAVE  TEMPERATURES  FOR  STORAGE  DIFFERENCES 


3085  S ( 1 0 ) = D ( 1 2 ) \S < 1 1 ) ==D(  9 ) \S ( 1 2 ) = D ( 1 0 ) \REM  TEMPERATURES 
3090  S(14)-^Ii(10)-D(9) 

3140  FOR  1=1  TO  14\T(I)=I(I)+S(I)\REM  SUM  FOR  DAY.  PRINT  HR  DATA 
3150  NEXT 

3160  S8-S(10)\S9=S(12)\RETURN  \REH  SAVE  FOR  TEMP  DIFFERENCE  - TUBE^HOUSE 
3200REM  + + + CALC  AVERAGES  ?.  FINAL  QUANTITIES 

3210  FOR  1=1 OTO  14\  REM  AVERAGE  TEMPERATURES  (NOTE  HEATS  NOT  INCL.,) 

3220  T(I)=T(I)/NO\REM  DAILY  AVGS  = SUM  / NR  HRS 
3230  NEXTXRETURN 

3400REM  SUBR  +++  PRINT  REPORT  HEADING 

ERIDAN  HOUSE  ( ID=  ' " ,R$ , " ^ " 


DELTA  ELECT  FURNC  HOUSE", 


3405 

'!f1\'«1"  DAILY  PERFORMANC 

E SUMMARY  FOR 

3406 

•Hn'WI"  UNITS:  ENERGY-MJ 

.TEMPERATURE- 

3410 

'»1\'tt1 ."REC« 

MO  DAY 

HR  " 

.TAB(16) 

f 

3420 

'«1."  SOLAR 

COLL 

STORE 

STORE 

DHUC 

3430 

'U1 ."  STORE 

AMBNT 

HOUSE 

STORE 

DEGRE 

3440 

!«1 ,TAB(16), 

3450 

'«1,"  INSO 

OUTPT 

OUTPT 

HOUSE 

OUTPT 

3460 

'111,  " TEMP 

TEMP 

TCMp 

LOSS 

HOUR 

3470 

'HI  ,TAB(16) , 

3495 

RETURN 

DHU 


POWER  DUCT 


LOAD 


Figure  8:  Computer  program  of  heat  flow  equations 
used  for  performance  summary 


29 


a: 

CD  cc 
uj  3 ^ 

o o •>- 

X 


TiJ 

<U 

•Vi 

n5 


ID 

a:  to 

O CO 


rO 


CO  _J  


tiO 

c 

*H 

0) 

rO 


UJ 

(n  0. 
rD  ar 
o UJ 
X H- 


I 


CO 


-p 

cj 

<D 


X 0. 
M X 
X UJ 

<C  I— 


II 

A 


UJ 

cn 


<c 


cfe: 

Ui 

X 

tn 

QC 

O 

u. 


o: 

<c 

X 

x: 

X 

to 


CJ 

CO 

UJ 

n 

I 

UJ 

Qd 

X 

H- 

<C 

QC 

UJ 

o. 

:c 

UJ 


UJ 

CJ  x 

X X 

<C  I 

X >- 

X CD 

o or 

U.  UJ 

or  X 

UI  UI 

o. 

>-  CO 


<C  X 

n X 


UI 

or  0. 
O X 

“o' 

03 

> 

*H 

• 

0 

rH 

H-  UJ 

pi 

CO  ►— 

N—- 

cd 

o 

tiO 

■r-3 

• 

UJ 

03 

a 

0 

CO  o 

• 

sd 

Pl 

'X  <X 

cr\. 

(U 

• 

bo 

•H 

o o 

rH 

0) 

CO 

0 

0 

X ^ 

P) 

rH 

•H 

£ 

rH 

o 

pi 

0 

•r-J 

o 

>» 

*• 

o 

o 

•o 

+» 

Pi 

cd 

•ri 

X 

0 

X o 

CO  . 

tOD 

(d 

■P 

0 

X X 

<0 

tlO 

• 

o 

0 

u.  n 

£ 

0) 

£ 

3 

0 

rH 

V 

& 

q 

►-3 

• 

0 

CQ 

•» 

o 

s 

0 

Td 

or 

Pi 

CO 

■I~3 

• 

• 

CJ  UI 

O 

•iH 

cd 

0 

OD 

•H 

0 

• 

UJ  X 

r- 

6 

+» 

ro 

X 

rH 

CO 

0 

q 

o 

Pi 

^ o 

M 

o 

+» 

Eh 

t»o 

• 

pi 

• 

rH 

•H 

II 

•iH 

UI  a. 

. 

Eh 

<0 

0 

0 

o 

0 

0 

a 

ft 

rH 

0) 

£ 

• 

rH 

•o 

o 

1 — 1 

s 

M 

rH 

& 

• 

0 

Pl 

a 

Td 

0 

0 

^cc 

« 

O 

(d 

£ 

«» 

rH 

O 

q 

0 

o 

II 

Td 

LC^ 

.J  X 

O 

w 

o 

Pi 

03 

0 

q 

•o 

bo 

a 

0 

, — ^ 

•H 

LPi 

UJ  o 

o 

-P 

in 

• 

O 

a 

0 

• 

0 

0 

CM 

0 

• 

X 

w 

(0 

-p 

CD 

q 

0 

■r-3 

E 

Ui 

Ui 

Pi 

0 

N. ✓ 

P> 

o 

ft 

Xi 

CO 

>5 

o 

rH 

cd 

bo 

Pi 

pi 

to 

0 

pl 

i— 

1— 

+=■ 

CO 

X 

pi 

0 

•> 

pi 

•H 

0 

Jh 

Pi 

O 

CJ  o. 

<D 

o 

t»D 

£ 

O 

0 

Td 

to 

0 

II 

X K- 

>> 

03 

0 

•l-J 

0 

0 

1 — f 

0 

X 

Td 

X X 

, 

X 

-p 

& 

X 

£ 

•» 

■P 

T 

0 

pj 

Td 

q 

iTd 

X CD 

V 

Cd 

■p 

cd 

0 

a 

0 

CJ 

-H 

a 

Eh 

o 

Pi 

to 

m 

£ 

£ 

0 

p: 

+» 

ft 

v 

-p 

O 

o 

0 

O 

•H 

u 

0 

Td 

•H 

P! 

0 

UJ  UJ 

-P 

-P 

H^ 

£ 

X 

X 

X 

s 

to 

0 

0 

0 

0 

Q 

CC  CO 
O X 
X 

ft 

(V 

Td 

0) 

CO 

0 

•b 

a 

-p 

o 

0 

Td 

0 

Pi 

Pi 

5 

Td 

0 

•H 

O 

pl 

o 

5 

O 

CO  X 

o 

Pi 

•p 

M 

M 

-p 

ft 

to 

0 

a 

h 

•H 

X 

•b 

Pi 

O) 

cd 

cd 

S’ 

•S 

ft 

0 

B 

Ch 

O 

0) 

> 

03 

Pi 

(3 

a 

Td 

0 

0 

Ch 

Pl 

UI  H- 

■P 

•H 

X 

o 

+> 

cd 

bo 

£ 

a 

0 

0 

X o. 

rH 

-p 

Pi 

-p 

•k 

0 

O 

0 

O K- 

m 

•H 

0) 

0) 

W 

% 

o 

sd 

0 

o 

0 

•b 

E 

O 

IS 

:3 

H-  X 

Td 

X 

+» 

•H 

o 

a 

to 

Pi 

-p 

CO  o 

c: 

nd 

-p 

£ 

ft 

Ui 

;s 

a 

a 

•H 

•b 

0 

0 

o 

o 

o 

p 

Pi 

Pi 

a 

Pi 

Ui 

X 

m 

'H 

C 

• 

Pi 

o 

§ 

rH 

o 

o 

•H 

O 

cn 

H- 

■P 

ro 

o 

0 

«w 

•p 

«H 

H 

-p 

+» 

a 

rH 

0 

o 

^ o* 

(d 

rH 

ft 

Ch 

Td 

0 

P! 

o 

• 

^ H- 
C9  D5 

CM 

1 

•H 

Td 

pd 

Td 

Td 

Pi 

0 

0 

0 

0 

PJ 

T- 

X o 

■rJ 

0) 

w 

O 

0 

0 

0 

£ 

0 

■p 

•H 

0 

to 

0 

(d 

+» 

u 

•f-3 

Pi 

Pi 

o 

IS 

O 

a 

a 

X 

pi 

a 

Pl 

II 

Pi 

o 

03 

cd 

0 

0 

o 

Pi 

X 

0 

o 

Pi 

0 

QC 

o 

> 

> 

ft 

X 

§ 

X 

0 

Ch 

IT) 

<C  O 

h 

rH 

03 

fjD 

•H 

•H 

Td 

> 

Ch 

Eh 

-J  CO 

03 

rH 

0 

rH 

rH 

Pi 

o 

-P 

a 

Ch 

Ch 

Ch 

a 

•H 

ft 

O X 

rH 

o 

PJ 

£ 

0 

0 

0 

•H 

a 

o 

O 

O 

O 

Td 

CO  »-♦ 

o 

o 

•H 

Td 

Td 

a 

Pi 

0 

rH 

§3 

1 

CD 

•k 

0 

-p 

X 

0 

0 

0 

0 

ckr 

03 

-p 

■P 

o 

-P 

Pi 

Pi 

Pi 

a 

rH 

tiD 

& 

cd 

0 

0 

a 

pi 

q 

q 

Pi 

B 

•• 

nr 

(d 

Pi 

q 

0 

0 

•H 

rH 

0 

-P 

-p 

-P 

o 

-p 

0 

+» 

03 

CD 

Pi 

X 

X 

0 

a 

X 

a 

a 

a 

+» 

a 

-P 

O 

C! 

X 

O 

0 

■H 

Pi 

Pi 

Pi 

0 

Pl 

•H 

<c 

■P 

0) 

o 

+> 

Pi 

H 

bo 

rH 

rH 

0 

0 

0 

0 

0 

g 

A 

'i 

ta 

cd 

cd 

p 

a 

•H 

Ui 

ft 

ft 

ft 

-p 

ft 

3 

0) 

03 

0 

rH 

rH 

m 

•p 

X 

pi 

£ 

E 

£ 

a 

£ 

O 

X 

XI 

X 

o 

o 

o 

o 

B 

o 

0 

0 

0 

0 

0 

g 

X 

Eh 

Eh 

Eh  -P 

CO 

CO 

O 

Eh 

<5 

M 

Eh 

Eh 

Eh 

w 

Eh 

o 

96 

0 

X 

yO— s. 

^ — s 

yi— V 

y<— N. 

-p 

UJ 

' 

T— 

CM 

ro 

LPv 

Md 

CO 

ON 

O 

r- 

CM 

ro 

o 

X 

— ’ 

T~ 

0) 

o 


§ 


ft 

o 

«H 


•H 

% 

<0 


0) 

rH 


Eh 


a^ 

0) 

u 

s. 

•H 


400 


31 


Figure  11;  Graph  of  5 Days  Data  in  April 
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Figure  12:  Graph  of  5 Days  Data  in 
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HOURLY  PERFORMANCE  DATA  FOR  SHERIDAN  SOLAR  SYSTEM 
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NOTICE 

This  report  was  prepared  as  an  account  of  work  sponsored  hy  the  Energy 
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through  the  Alternative  Renewable  Energy  Sources  Program.  Neither  the  State 
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the  accuracy,  completeness  or  usefulness  of  any  information,  apparatus,  pro- 
duct or  process  disclosed,  or  represents  that  its  use  would  not  infringe 
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ABSTRACT 

This  project  was  a site-built,  active-air  solar  system  retrofit  to  an 
older  residence  in  Hamilton,  Montana.  The  grant  application,  design  and 
execution  were  carried  out  by  the  District  #11  Human  Resources  Development 
Council.  The  design  goal  vras  to  produce  a *’lov;  cost"  system.  The  result- 
ing design  was  unconventional  in  several  respects:  the  rock  bin  incorpora- 
ted a baffle  which  reduced  its  storage  capacity  by  a factor  of  two  and  a 
damper  was  not  used  to  isolate  the  collector  loop  and  the  house  loop. 

These  design  anomalous  plus  numerous  air  leaks  reduced  the  potential  per- 
formance of  this  system.  The  system  supplied  negligible  solar  heat  to  the 
house  during  the  monitoring  period.  It  seems  that  the  system  sometimes 
removes  heat  from  the  house. 


NOTE  ON  UNITS 

The  data  is  presented  in  metric  units.  These  units  are  in  accord  with 
the  recommendations  of  the  Intemational  Solar  Energy  Society  as  well  as 
most  current  technical  journals.  The  overall  siunmaries  and  numbers  in  the 
test  have  been  given  in  dual  units.  The  heat  load  calculations  are  pre- 
sented in  conventional  units. 

Most  readers  are  familiar  with  the  Celsius  or  Centigrade  temperature 
scale  and  are  also  familiar  with  energy  units  of  kilowatt  hours.  The  only 
other  metric  unit  which  occurs  frequently  in  this  report  is  the  Joule  (j), 
the  metric  unit  of  energy  v/hich  replaces  the  Btu  (a  Btu  is  equal  to  1,055 

Joules).  The  mega  joule  (MJ)  is  equal  to  lO^J  (100,000  Btu/hr  = IO5  Mj/hr). 

2 

Other  abbreviations  used  in  this  report  include:  m = meters;  m = 
square  meters;  1 = liters;  h = hours;  gph  = gallons  per  hour;  and  Iiffitu  = 
millions  of  Btu’s. 
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1.0  INTRODUCTION 

This  project  is  located  in  Hamilton,  Montana  and  is  the  residence  of 
Mrs.  Pauline  Salyer  and  her  sons  Jay  and  Robin.  This  house  was  chosen  by 
the  District  #11  Hiiman  Resources  Development  Coiincil  to  demonstrate  retro- 
fitting an  older  house  with  a solar  heating  system.  The  goal  of  their  pro- 
gram was  to  "cut  use  of  conventional  heating  energy  by  at  least  50^  "the 
target  homes"  and,  thus,  reduce  the  economic  burden  of  rising  energy  costs 
on  low-income  persons. 

The  proposed  project  was  submitted  by  Ken  Boggs  and  Roger  King  and  six 
solar  systems  were  installed  on  low-income  homes.  Labor  to  build  the  sys- 
tem was  supplemented  with  Comprehensive  Employment  and  Training  Act  (CETA) 
workers. 

2.0  ACKNOWLEDGMENTS 

We  gratefully  acknowledge  the  cooperation  of  Mrs.  Pauline  Salyer  during 
this  monitoring  project.  Mrs.  Salyer  allowed  us  to  use  one  of  her  closets  to 
house  the  data  acquisition  system.  V/e  spent  several  days  crawling  around  in 
the  attic  and  crawl  space  measuring  flow  rates  and  locating  temperature 
transducers.  During  the  installation,  dozens  of  trips  were  made  in  and  out 
of  the  house  and  her  kitchen  and  bedroom  were  often  littered  with  instru- 
ments, probes  and  meters.  Mrs.  Salyer  graciously  accepted  these  disruptions 
and  was  always  helpful. 

Roger  King  is  acknowledged  for  spending  several  hours  describing  the 
construction  of  the  solar  system. 

3.0  PROJECT  DESCRIPTION 

Photographs  of  the  Salyer  house  and  the  solar  system  are  shown  in  Fig- 
ures 1 , 2 and  3.  The  floor  plan  of  the  house  and  a cutaway  of  the  air 
collector  construction  is  shown  in  Figure  4 and  a schematic  of  the  solar 
system  is  shown  in  Figure  5» 

2 

The  collector  is  air-cooled  and  has  a total  aperture  of  28.46  m (309 

2 

ft  ) . The  collector  is  single  glazed  with  tempered  patio-door  panels  0.3  cm 
(0.125  in)  thick.  The  absorber  is  thin  aluminxim  painted  flat  black  and  the 
panels  are  insulated  with  plastic  foam.  Figure  3.  Air  flows  behind  the  ab- 
sorber and  baffles  are  used  to  distribute  the  flow.  The  panels  were  site- 
built  and  permanently  attached  to  the  roof.  The  roof  faces  south  and  is 
pitched  at  a 45^  angle. 
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The  schematic,  Figure  5i  shows  a two-fan  system  with  a single  motorized 
damper  in  the  collector  loop.  A differential  thermostat  was  installed  to 
control  the  collector  fan  and  damper.  This  control  has  been  disabled  and, 
during  the  monitoring  period,  the  collector  fan  and  damper  were  controlled 
by  a simple  thermostat  located  behind  the  solar  absorber.  This  thermostat 
was  set  at  100°P.  The  house  fan  is  controlled  by  a thermostat  within  the 
living  space,  A second  thermostat  in  the  rock  bin  heat  storage  opens  the 
house  fan  circuit  if  the  rock  bin  is  cold,  thus  preventing  cold  storage  air 
from  entering  the  house, 

A roof  vent  with  a manual  damper  is  intended  for  use  during  the  summer 
to  vent  excess  collector  heat.  Figures  2 and  5»  The  ductwork  consists  of 
insulated,  circular  flexible  duct  combined  with  sections  fabricated  with 
2.5  cm  (1  in)  thick  fiberglass  duct  board  (see  photographs  in  Figure  3) • 

The  designer  chose  not  to  use  a manufactured  solar  air  handling  unit 
to  reduce  the  cost  of  the  system.  A common  duct  above  the  ceiling.  Figures 
3 and  5i  is  intended  to  be  a cold  air  supply  to  the  collectors  as  X’jell  as  a 
cold  air  return  from  the  house.  There  is  no  damper  to  positively  close  off 
the  house  loop  which  may  explain  the  unconventional  baffle  in  the  top  of 
the  rock  bin.  Figure  5» 

Axixiliary  heat  is  supplied  to  the  house  by  a circulating  natural  gas 
heater  controlled  by  an  integral  thermostat.  This  heater  is  located  in 
the  dining  room.  Figure  4»  There  is  also  a x^rood  stove  on  the  enclosed  porch 
next  to  the  kitchen  which  is  used  for  supplementary  heat, 

4.0  INSTRmENTATION  AND  MEASUREtffiNTS 

4. 1 Transducer  Layout 

The  transducer  layout  is  shown  in  Figure  Solar  insolation  was  mea- 
sured in  the  plane  of  the  collectors  (45°)  at  mid-point  of  the  western  edge 
of  the  array.  The  ambient  temperature  was  measured  by  a shaded  transducer 
on  the  north  side  of  the  house  and  about  3 m (IO  ft)  above  the  groxind,  Tvjo 
temperature  probes  were  mounted  inside  the  house  and  one  was  placed  in  the 
attic. 

The  solar  collectors  had  two  independent  probes  in  both  the  inlet  and 
outlet  ducts  near  the  collectors.  Two  probes  were  used  to  provide  backup 
as  well  as  a check  point  on  this  important  measurement.  The  rock  bin— to— 
house  loop  had  inlet  and  outlet  duct  probes  to  measure  stored  solar  heat  en- 
tering the  house.  A status  relay  was  connected  to  each  fan  to  record  when 
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the  fan  was  running.  The  computer  program  contained  logic  to  zero  the  duct 
temperatures  if  the  fan  was  off.  This  allows  calculation  of  the  true  duct 
temperature  while  the  fan  is  running  in  a cyclic  mode. 

A temperature  probe  was  inserted  in  the  hot  end  of  the  rock  bin  storaige. 
Another  probe  was  suspended  in  the  insulated,  unheated  shed  which  contains 
the  rock  bin.  The  total  electric  power  was  monitored  with  clamp-on  ammeters. 
These  readings  were  calibrated  against  the  service  meter. 

The  heat  output  of  the  free  standing  wood  stove  and  axixilicLry  natural 
gas  heater  could  not  be  accurately  monitored.  The  normal  way  to  monitor 
this  type  of  heating  device  is  to  enclose  it  in  a calorimeter  room.  This 
can  be  done  in  a laboratoiy,  but  not  in  an  occupied  residence.  Readings  of 
the  natural  gas  meter  during  the  test  period  allowed  calculation  of  the 
average  auxiliary  heat. 

4.2  Plow  Measurements 

Figure  6 shows  some  typical  results  of  measurements  of  the  air  flow  in 
the  ducts.  A hot-HfJire  anemometer  was  used  to  measure  the  mass-flow  velocity 
at  a network  of  points  on  the  cross  section  of  each  duct.  These  results 
were  then  integrated  to  yield  net  flov:  through  the  duct. 

To  check  the  overall  accuracy  of  the  technique,  the  flow  was  measured 
at  the  collector  fan  inlet  and  outlet  for  one  mode.  Figure  6 shows  that 
these  results  compare  within  This  is  considered  as  satisfactory  agree- 

ment considering  the  turbulence  and  short  lengths  of  duct  at  the  measurement 
stations. 

4»  3 Flow  Balance  and  Leaks 

Probably  the  most  significant  result  of  this  monitoring  effort  is 
shown  in  Figures  7 and  8.  In  Figure  7 the  collector  fan  was  on  and  the 
house  fan  was  off.  The  numbers  by  the  arrows  are  volumetric  flow  rates  in 
cubic  meters  per  minute  (m^/min).  The  nxxmbers  represent  measurements  except 
when  enclosed  in  parentheses.  The  enclosed  numbers  are  leakage  flows  neces- 
sary to  balance  the  measured  flow  rates. 

Looking  first  at  the  rock  bin  heat  storage  box  we  see  that  13.2  m^/min 
axe  flowing  into  the  top  of  the  box.  Ideally,  this  same  flow  should  be 
leaving  the  bottom  of  the  box  and  proceeding  back  up  the  cold  air  return 
duct  to  the  collectors.  Measurements  show  that  only  0.8  m^/min  is  returning 
through  the  cold  air  duct  and  that  3.7  m^/min  is  flowing  through  the  house 
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duct.  Flow  balance  of  the  inputs  and  outputs  of  the  rock  storage  imply  that 
8.6  m^/min  is  leaking  out  of  the  storage  bin. 

Looking  at  the  collector  array  in  Figure  7i  we  measured  a flow  of  2.2 
m^/min  leaking  through  the  summer  vent  (which  should  be  closed) . The  flow 
in  the  return  vents  in  the  ceiling  of  the  living  room  and  bedroom  was  mea- 
sured as  approximately  4 to  6 m^/min.  These  flows  should  also  be  zero. 

Flow  balance  on  the  collector  implies  that  leaks  amounting  to  about  4*2  m^/ 
min  exist  in  the  collectors  and/or  their  ducting. 

Figure  8 shows  results  of  another  set  of  measurements  of  the  flows  in 
the  duct  system  with  both  collector  and  house  fans  on.  These  measurements 
were  more  approximate  than  those  of  Figure  7»  The  main  conclusions  of  this 
measurement  set  are  that  (l)  large  rock  bin  and  house/collector  leakages 
exist,  and  (2)  there  is  almost  no  flov;  in  the  collector  cold  air  return. 

It  appears  that  the  pressure  drop  through  the  thickness  of  the  rock 
bin  is  much  greater  than  the  pressure  drop  across  the  internal  baffle.  The 
air  that  enters  the  rock  bin  either  leaks  out  or  flov/s  around  the  baffle  and 
into  the  house. 

4. 4 Data  Reduction 

The  existence  of  the  unexpected  and  unusual  flow  patterns  and  large 
leaks  in  this  system  makes  precise  analysis  of  the  data  impossible.  V.'hen 
the  collector  is  pres\imably  storing  heat  in  the  rock  bin  and  there  is  pre- 
sumably no  flow  throu^  the  house  we  find  that  collector  air  that  is  not 
lost  due  to  leakage  is  actually  short-circuiting  through  the  house.  The 
collector  itself  which  should  be  drawing  air  from  the  rock  bin  is  actually 
drawing  in  air  from  (a)  the  house,  (b)  the  summer  vent,  and  (c)  some  place 
else,  probably  distributed  leaks. 

Collector  efficiency  is  defined  in  terms  of  a device  which  has  one  in- 
let and  one  outlet.  This  collector  seems  to  be  drawing  air  from  at  least 
four  separate  sources. 

As  a ’’best  effort”  attempt  to  analyze  the  data,  a weighted  average 
collector  inlet  temperature  was  created  based  on  the  flow  pattern  in  Figure 
7;  165^  was  assumed  to  be  at  ambient  temperature,  32^  at  attic  temperature, 
45^  at  average  house  temperature  and  7^  at  rock  bin  return  temperature.  The 
temperature  of  the  leakage  collector  air  entering  the  house  was  assumed  to 
be  3°C  below  the  air  temperature  entering  the  storage.  This  estimate  was 
necessary  because  the  data  system  ignored  the  house  duct  temperature  because 
the  house  fan  was  never  on. 
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During  the  test  period  the  house  fan  never  did  turn  on,  due  probably 
to  low  storage  temperatures.  Thus,  the  only  solar  heat  that  entered  the 
house  was  the  estimated  heat  due  to  flow  leaking  around  the  internal  baffle 
in  the  rock  bin, 

4, 5 Data  Base 

The  data  acquisition  system  was  installed  on-site  and  operated  from 
January  21,  1979  to  March  l8,  1979*  Eighteen  channels  of  data  were  re- 
corded, Unfortunately,  the  second  cassette  tape  containing  data  from 
February  12  to  March  18  could  not  be  read.  It  had  evidently  been  exposed 
to  a magnetic  field  of  unknown  origin. 

The  data  base,  thus,  consists  of  9180  data  points  taken  during  510 
hours  or  21  days  betv;een  January  21  and  February  11,  1979-  Each  data  point 
represents  about  600  averaged  readings,  thus,  there  were  about  5i 500,000 
individual  readings, 

5,0  HOURLY  MD  DAILY  SUMMAEY  OF  DATA 

Figure  9 shov/s  the  equations  used  to  format  the  summary  data  for  this 
project.  The  hourly  performance  summary  is  shown  in  the  tables  in  Appendix 
I,  Table  1 is  included  in  the  main  text  and  will  be  discussed  as  an  example. 
Units  of  energy  are  mega  joules,  MJ,  and  temperatures  are  in  degrees  Celsius, 
°C,  Table  1 shows  the  behavior  of  the  system  on  a typical  cold  and  sunny 
day.  The  sun  strikes  the  collector  at  9 AM  but  there  are  apparently  some 
clouds  until  10  AM,  Between  1 1 AM  and  noon  the  collector  fan  turns  on  and 
stays  on  until  shortly  after  3 PM, 

During  this  time  the  collector  outputs  approximately  J2  MJ  (76,000  BTU) 
of  which  8 MJ  (8,440  BTU)  leaks  into  the  house.  The  hourly  heat  load  of  the 
house  is  around  25  Mj/hr  (26,400  BTU/hr)  so  the  solar  heat  gain  is  veiy  small. 
Note  that  the  temperature  of  the  shed  i^hich  contains  the  rock  bin  increases 
2°c/hr  (3.6°F/hr)  while  the  collector  fan  is  on.  This  is  probably  due  to  the 
hot  air  leak  discussed  previously. 

The  storage  bin  temperature,  while  the  collectors  are  on,  reflects  ba- 
sically the  outlet  air  temperature  of  the  solar  collectors.  This  reaches  a 
maximum  of  40,5°C  (105°F)  which  is  about  typical  for  this  system  over  the 
cold  test  period.  This  performance  is  not  remarkable  since  the  collectors 
are  sucking  much  of  their  inlet  air  out  of  the  house  vjhich  is  only  about  1 5°C 
(27°F)  below  their  outlet  temperature,  VIhen  the  fan  turns  off,  the  tempera- 
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ture  at  the  top  of  the  rock  bin  is  30 °C  (86^)  and  drops  to  25°C  (77°F)  by 
midnight.  This  temperature  is  too  low  to  use  in  the  house  and  the  energy 
is  dissipated  during  the  night. 

The  attic  temperature  shows  an  increase  of  several  degrees  on  some 
mornings  ciround  7 AM.  This  may  be  due  to  a heat  leak  from  the  auxiliary 
heater.  The  bottom  line  in  the  table  gives  the  sum  of  the  hourly  energy 
flows  and  the  average  temperatures. 

6.0  OVERALL  PERFORIiANCE  SMAKY 

Table  2 summarizes  the  hourly  data  given  in  Appendix  I,  The  hourly 
heat  quantities  are  summed  over  the  entire  test  period  and  temperatures  are 
averaged.  The  first  part  of  Table  2 shows  that  2^  of  the  solar  energy 
striking  the  collector  surfaces  was  delivered  to  the  house.  This  heat  was 
entirely  due  to  the  leakage  flow  around  the  baffle  in  the  rock  storage.  The 
house  fan  did  not  turn  on  during  the  test  period  so  no  stored  solar  heat  was 
delivered  to  the  house. 

The  energy  inputs  are  compared  to  the  calculated  heat  load  in  the  sec- 
ond part  of  Table  2,  The  solar  system  is  seen  to  contribute  a negligible 
heat  input,  amounting  to  only  0,6^  of  the  load.  The  electrical  energy  to 
run  the  fan  used  about  l/3  of  this  amount  of  energy  and  was  dissipated  in 
the  attic. 

The  house  temperature  averaged  25°C  (77°F)  during  the  test  period  while 
the  ambient  temperature  averaged  -7.4°C  (l8.7*^P)«  The  temperature  of  the 
hot  (top)  end  of  the  rock  bin  averaged  22,4°C  (72°P)  which  is  clearly  belov; 
the  house  temperature  and  not  usable. 

7.0  COLLECTOR  EFFICIENCY  AND  SHADING 

The  collector  efficiency  curves  for  this  system  are  shown  in  Figure  10. 
The  collector  output  as  well  as  the  inlet  temperature  v;ere  calculated  using 
the  same  assiimptions  mentioned  previously.  Based  on  this  approach,  the  ef- 
ficiency falls  v;ell  belov:  the  HUI)  reference  for  this  type  of  collector. 

This  seems  plausible  as  the  collector  inlet  temperature  is  high  ( since  it 
draws  heated  air  from  the  house)  and  since  there  are  sizable  leaks  from  the 
vent  and  the  attic  which  dilute  the  solar-heated  air. 

Figure  11  shov:s  the  shading  of  the  Salyer  collector  during  the  winter. 
There  is  some  early  morning  shading  due  to  trees  along  the  street  in  front 
of  the  house  and  afternoon  shading  due  to  trees  and  mountains  to  the  west. 
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The  collector  is  exposed  to  the  sxm  during  the  middle  of  the  day  when  the 
radiation  is  most  intense  and  usable.  It  can  be  concluded  that  these  col- 
lectors have  good  exposure  to  the  sun  and  that  the  shading  would  cause  in- 
significant losses  in  efficiency. 

8.0  HEAT  LOAD  CALCULATIONS 

The  calculations  of  heat  load  for  the  Salyer  house  are  shown  in  Table 
3.  Part  of  the  ceiling  of  this  house  is  covered  vjith  an  attic  having  ver- 
miculite  and  fiberglass  batts  producing  an  insulation  value  of  R-38,  The 
roof  over  the  kitchen  and  the  back  porch  is  a shed  roof  v/ith  R-1 1 insulation. 
Most  of  the  windows  have  storm  windows.  The  largest  item  in  the  heat  load 
is  the  walls  which  are  not  insulated.  The  heat  load  factor  used  in  the  per- 
formance summary  was  reduced  from  the  calculated  value  by  25^.  This  was 
done  since  the  bedrooms  and  back  porch  are  often  closed  off  and  are  remote 
from  the  auxiliary  heater  in  the  living  room.  This  produced  a value  of  heat 
load  which  fit  the  experimental  data  better, 

9.0  PROJECTED  POTENTIAL  PERFORMANCE 

Figure  12  contains  a plot  of  the  output  of  the  solar  collector  for  a 
variety  of  values  of  daily  total  insolation.  These  values  were  taken  from 
the  siimmary  data  in  Appendix  I,  This  characteristic  is  combined  with  1978- 
79  solar  radiation  data  from  the  Hamilton  High  School  and  long-term  degree 
day  for  Hamilton  to  produce  the  idealized  annual  performance  projection  shown 
in  Table  4.  This  projection  assumes  that  the  present  system  is  modified  so 
that  the  heat  delivered  to  the  rock  bin  is  actually  used  in  the  house.  This 

analysis  predicts  an  annual  solar  fraction  of  22^  and  an  overall  collector 
efficiency  of  13^» 

The  last  column  in  Table  4 is  auxiliary  heat  calculated  from  Mrs,  Sal- 
yer *s  utility  bills.  Through  November  these  values  compare  reasonably  well 
with  calculated  heat  loads  using  average  degree  days.  During  the  winter,  the 
fuel  usage  is  over  60^  higher  than  average. 

Mrs.  Salyer  reports  that  last  winter  her  fuel  bills  were  a lot  higher 
than  normal.  I believe  she  asked  the  utility  company  to  replace  her  meter 
because  she  thought  it  v/as  reading  high.  There  are  three  factors  which  could 
contribute  to  these  higher  bills:  (a)  last  winter  was  exceptionally  cold, 

(b)  there  were  utility  rate  increases,  and  (c)  the  solar  system  draws  hot  air 
from  the  house  and  deposits  significant  amounts  of  heat  into  storage  where  it 
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is  lost.  There  is  not  enough  data  to  determine  the  relative  values  of  these 
factors  at  this  time. 

10.0  GRAPHICAL  PRESENTATION  OF  DATA 

Graphs  of  ten  da;ys  of  data  are  shown  in  Figures  13  and  14.  The  ambient 
temperatures  were  very  low  during  the  first  seven  days  of  this  period.  Six 
of  these  days  were  clear  and  the  collectors  operated  for  several  hours  each 
day.  The  daily  peaks  in  the  temperature  of  the  shed  containing  the  rock-bin 
heat  storage  are  due  to  heat  loss  in  this  device.  The  attic  temperature  is 
approximately  midway  between  house  and  ambient  temperature.  Since  the  attic- 
house  interface  is  insulated  to  R-38,  this  high  attic  temperature  indicates 
that  hot  air  from  the  house  is  entering  the  attic.  The  house  temperature 
shows  a small  decrease  each  day  when  the  solar  collector  fan  turns  on.  This 
behavior  is  consistent  with  the  flow  measurements  depicted  in  Figure  7 which 
shows  that  the  collector  fan  pumps  air  from  the  house. 

11.0  CONCLUSIONS 

This  solar  system  did  not  produce  significant  heat  during  the  cold 
test  period.  If  no  modifications  are  made,  the  system  should  be  turned  off 
during  cold  weather.  Modifications  which  appear  advisable  to  increase  the 
energy  efficiency  of  this  house  are  listed  below  in  the  order  of  importance: 

(1)  Add  motorized  damper  to  house  loop  and  remove  baffle  from  heat 
storage  bin. 

(2)  Seal  air  leaks. 

(3)  Lower  the  house  temperature. 

(4)  Close  summer  vent  securely, 

( 5)  Increase  insulation  of  heat  storage  bin. 

(6)  Seal  off  ceiling  vents  and  add  floor— level  cold  air  return  to 
rock  bin, 

(?)  Adjust  and  use  the  differential  thermostat  and  check  operation  of 
other  controls, 

(8)  Equalize  flow  in  collectors. 


TABLE  1 

EXAJIPLE  OF  HOURLY  PERFORMANCE  DATA 
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TABLE  2 

OVERALL  PERFORMANCE  SUMI4ARY  OP  SALYER  HOUSE  DURING  MONITORING  PERIOD 


^IBTU 

Total  Insolation 

6750 

7.1 

100 

Collector  Output 

986 

1 

15 

Solar  to  House 

96 

0.1 

2 

Solar  Storage  to  House 

0 

0 

0 

*House  Heat  Load 

(14142) 

14.9 

100 

5 

‘Solar  Heat 

96 

0.1 

6 

Electric  Input 

3315 

3.5 

23 

^Aixxiliary 

(10731) 

11.3 

76 

Average  Temperatures 

House 

25 

77 

Ambient 

- 7.4 

18.7 

Rock  Bin  Storage  (Top) 

22.4 

72 

Calculated,  -20^^ 

2 

Note:  The  fan  absorbed  about  35  MJ  of  electricity  during  the  test  period, 

^ Mega  Joules 
4 


Millions  of  BTU*s 
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TABLE  3 

CALCULATED  HEAT  LOAD  FACTOR  FOR  SALYER  HOUSE 


R 

U 

AREA 

U X A 

(Btu/hr  ft^  °F) 

(sq.  ft.) 

Roof 

Shed 

11  .091 

360 

33 

Attic 

38  .026 

754 

20 

V/alls 

Nominal 

4 .25 

840 

210 

Rock  Bin 

24  .041 

102 

4 

Floor 

^8  .12 

1114 

133 

VI  indows 

Single 

.9  1.13 

24 

27 

Storm 

1.8  .55 

186 

103 

Infiltration: 

10,000  X 0.5  X 0.018 

90 

Total  Heat  Load  619  Btu/hr°F 

or 

1.17  Mj/hr°C 

Effective  R,  over  crawl  space 
Assumin^j  -g-  air  change  per  hour 


PROJECTED  POTENTIAL  YEARLY  PERFORMANCE 
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This  assumes  system  is  modified  to  increase  efficiency 
Based  on  utility  bills. 


(a)  East  side  of  Salyer  house 


V 


(h)  North  side  of  Salyer  roof  showing 
collectors  extending  past  ridge. 


Figure  1:  Salyer  house  photographs 
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(a)  Solar  collectors,  vent  and  rock 
bin  enclosure 


(c)  Close  up  of  collectors 


(b)  Plywood  rock  bin  inside 


attached  enclosure 


( d)  SE  corner  of  house 





Ti'i  omT*o  O • 
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(a)  Collector  supply  and  return  ducting 


(b)  Collector  fan,  ducting  and  wiring 


Figure  3:  Attic  of  Salyer  house 


Flashing 


^ l/8"  Tempered  Glass 


Absorber 


Baffle 

(polyi  so-synur^ite ) 
nermax  Insulatior' 


f 

\ 


\ 

.020"  thick  (Flat  Black) 


^Existing  Roof 


Figure  4s  Floor  plan  of  house  and  schematic  of  collector  construction 
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Ambient 


X 

/' 


Collector ; 
Fan  : 


i • 


Motorized'  | 
Damper  ! 


Baffle 


Rock 

Bin 

Heat 

Storage 


I • • 


House 


Fan 


Collectors 


-._j 


House 


Power  Input 


"Living  room 


•Attic 


• Temperature  Probe 


^ Status 


Figure  5s  Schematic  of  solar  heating  ^stem 

showing  instrumentation  arrangement. 
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30.5  cm 


\ 


r- 

0 


1 2 

House  Duct,  House  Faji  On  : Flow=  13.54  m^/min,  T=l8° 


k 19cm  > 

i ■ 

46 

cm 

r_.. 

\ 

\ 

\ 


Fan  Outlet 


Collector  Out,  Collector  Fan  On  : Flow  = 16.35  m^/min,  T=28°C 


6 m/e 


^ tn/s 


\ 


30. 5cm 


Pan  Inlet 

Collector  Out,  Collector  Pan  On 


Plow-  15.24  mVrain,  T-280 


6 m/e 


Figure  6:  Plow  profiles  resulting  from  scanning  the 
ducts  using  a hot  wire  anemometer. 
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Figure  7:  Flow  Balance  and  LeaJcs;  Collector  fan  - on, 
House  fan  - off.  Flow  in  m’^/min. 
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Collector  Arra^y 


House 


Figure  8:  Plow  Balance  and  Leaks j Collector  fan  is  on  and 
House  fan  is  on*  Flow  in  cubic  meters/minute* 
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JOOOREM  SUER  + »+  COMPUTE  HOURLY  SUilS  + + t + f t+.+ + 

3005  "'--  SU««IHG" 

3010  iVO  = HO  + l\RE«  HHRS 

3012  RE«  ALL  THE  FOLLOLIING  ARE  I«  MEGA  JOULES 
3015  S(1)  = <U(1  )-.01  )t.«t)  + 2B.46  + 3.6\REf1  SOLAR”lNPUr 

3018  T1  = .U  + D(7)  + .32*Ei(10)  + . 45+ ( » ( 8 ) + [!  ( 9))/2  + .07*D(11)\REM  AUG  TIN 
3020  S(2)  = 13.2+.06  + 1 .2  + 1 .01  + ((ri(15)/(D(4)  + .0001  ))-T1)+ti(4) 

3025  T2=(Df8)+3+D(9) )/4  \ REN  AVERAGE  HOUSE  TEMPERATURE 
3030  S(3)=3.7+.0A+1 .2+1 .01*( (D(15)/(D(4)+.0001 ))~3-T2)  +D(4) 

3035  REM  ABOVE  IS  SOLAR  HEAT  TO  HOUSE.  ASSUMES  3C  TEMP  DROP 

3040  S(4)=Ii(5)  + 16.A+.06  + 1 .2  + 1 + ([!(  9 ) -D  ( 1 3 ) ) \REM  STORAGE  HEAT  TO  HOUSE 

3050  S(5)  = ([i(2)+D(3)  ) + 3.6  + . 494  \REM  ELECTRICAL  POUER  INPUT 

3060  S(6)=1 .17+(T2-D(7))+.75  \ REM  HOUSE  HEAT  LOAD,  CORRECTED  X.75 

3062  S(7)=T2-D(7)  \ REM  DEGREE  HOURS 

3065  S(8)=(D(2)+D(3) )+.494  \REH  KU  INPUT,  CORRECTED 

3070  S ( 9 ) =D (1  2 ) ' S ( 1 0 ) =[!  n 6 ) \S < 1 1 ) = D (1 0 ) \S  ( 1 2 ) =T2\S ( 1 3 ) =D  ( 7 ) NREM  TEMPERATURES 
3080  S(1 1 )=ri(2)+D(3) 

3085  S(14)  = (D(n-.01  ) + .95+10 

3130  'B1  ,X3I,R,X4I,M0,3:3I.D0,X5I,H0,TAB(15) , 

3140  FOR  1 = 1 TO  14\T<n  = T(I)+S(I)\REM  SUM  EOR  DAY,  PRINT  HR  DATA 
3150  'til  .%7F1  ,S(I).\NEXT\'}f1 
3160  RETURNN'RI 

3200REM  +++  CALC  AVERAGES  X FINAL  QUANTITIES 

3210  FOR  1=  9T0  14\  REM  AVERAGE  TEMPERATURES  (NOTE  HEATS  NOT  INCL.,) 

3220  T(I)=T(I)/N0\REM  DAILY  AVGS  = SUM  / NR  HRS 
3230  NEXTNRETURN 

3400REM  SUBR  +++  PRINT  REPORT  HEADING 
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Figure  9:  Computer  equations  for  heat  quantities 
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Figure  10:  Collector  efficiency 
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Figure  11:  Winter  shading  of  Salyer  collector 
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Figure  12:  Daily  avereige  collector  performance 
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Figure  13;  Graph  of  5 Days  Data  in  January 
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Figure  14:  Graph  of  5 Da^ys  Data  in  Pehruaiy 
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ABSTRACT 

This  residence  is  ovmed  hy  the  George  Mattson  family  and  is  located 

about  three  miles  south  of  Bozeman,  Montana,  The  active  air  solar  system 
2 

uses  600  ft  of  collectors,  a single  fan  and  a rock  bin  for  heat  storage. 
The  performance  of  the  house  was  monitored  during  the  late  winter  and 
spring  of  1979*  This  is  a new  house  and  the  Mattsons  moved  into  it  in 
mid-April  1979* 

The  monitoring  data  showed  that  this  system  had  several  operational 
problems  including  air  leaks,  damper  malfunctions  and  erratic  control  be- 
havior, The  blov/er  motor  was  replaced  at  one  point.  The  owner  and  con- 
tractor are  in  the  process  of  making  major  changes  in  the  system  in  hopes 
of  improving  its  performance.  The  low  performance  of  this  project  reflects 
the  net  effect  of  these  problems. 

During  the  monitoring  period,  the  collectors  delivered  a total  of  15?^ 
of  the  available  solar  insolation  to  the  house.  Of  the  solar  heat  sent  to 
storage,  80^  was  lost.  The  annual  solar  fraction  was  estimated  to  be  a 
maximum  of  yyfo.  If  the  system  problems  can  be  fixed,  the  solar  efficiency 
should  be  about  twice  its  current  value. 


NOTE  ON  UNITS 

The  data  is  presented  in  metric  units.  These  xinits  are  in  accord  vzith 
the  recommendations  of  the  International  Solar  Energy  Society  as  well  as 
most  current  technical  journals.  The  overall  summaries  and  n\xmbers  in  the 
test  have  been  given  in  dual  units.  The  heat  load  calculations  are  pre- 
sented in  conventional  units. 

Most  readers  are  familiar  with  the  Celsius  or  Centigrade  temperature 
scale  and  are  also  familiar  with  energy  units  of  kilowatt  hours.  The  only 
other  metric  unit  which  occurs  frecjuently  in  this  report  is  the  Joule  (j), 
the  metric  unit  of  energy  which  replaces  the  Btu  (a  Btu  is  equal  to  1,055 

Joules).  The  mega  joule  (MJ)  is  equal  to  lO^J  (100,000  Btu/hr  = 105  Mj/hr). 

2 

Other  abbreviations  used  in  this  report  include:  m = meters;  m = 
square  meters;  1 = liters;  h = hours;  gph  = gallons  per  hour;  and  MBtu  = 
millions  of  Btu*s. 
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2.0  INTRODUCTION 

The  Mattson  residence  is  located  about  three  miles  south  of  Bozeman. 

The  house  and  the  solar  system  were  designed  by  Mr.  Mattson  who  is  a regis- 
tered architect.  The  Mattson  family  did  a large  portion  of  the  construction 
of  this  house  which  has  spanned  1977  to  the  present.  The  house  has  been  oc- 
cupied since  mid-April,  1979» 

3.0  DESCRIPTION  OF  HOUSE 

Figures  1 and  2 show  photographs  of  the  Mattson  solar  house  and  a floor 
plan  and  elevation  view  are  shown  in  Figures  3 and  4.  The  house  is  well  in- 
sulated and  is  oriented  to  conserve  energy.  Solar  heat  is  supplied  to  the 
house  by  an  active  air  system  combined  with  passive  clerestory  windows.  An 
internal  masonry  wall  receives  radiation  from  the  winter  sun  through  the 
clerestory  windows.  A natural  gas  furnace  and  a wood  stove  supply  auxiliary 
heating.  The  calculated  heat  load  for  this  house  is  shown  in  Table  1. 

A shading  diagram  for  a point  1.5  m (5  ft)  above  the  groinid  at  the  cen- 
ter of  the  collector  array  is  shown  in  Figure  5*  A large  tree  to  the  east 
of  the  house  is  seen  to  cause  some  shading  early  in  the  morning.  The  outline 
of  several  small  trees  along  the  property  boundary  to  the  south  are  shown  but 
do  not  cause  significant  shading  at  this  time.  These  trees  were  planted  by 
Mr.  Mattson *s  neighbors  after  the  house  was  built  in  order  to  block  objec- 
tionable reflections  from  the  collector. 
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4.0  THE  SOLAR  SYSTEM 

A schematic  and  control  matrix  of  the  Mattson  solar  system  is  shown  in 
Figure  6,  This  figure  also  shows  the  locations  of  transducers  used  to  moni- 
tor the  performance  of  this  house.  The  control  matrix  shows  the  positions 
of  the  three  dampers  for  the  primary  modes  of  operation. 

A single  fan  powers  the  system  and  three  motorized  dampers  configure  the 
flow.  A differential  thermostat  signals  the  availability  of  heat  from  the 
collectors.  A two-stage  thermostat  signals  demand  for  heat  by  the  house;  the 
first  stage  calls  for  solar  heat  and  the  second  stage  starts  the  auxiliary 
gas  furnace.  Hot  air  is  distributed  to  the  house  via  an  insulated  plenum  in 
the  crawl  space  and  vents  cut  through  the  floor.  All  ducts  are  insulated. 

The  collectors  are  double  glazed  and  are  integrated  into  the  structure 
of  the  south  wall  of  the  mechanical  room.  The  absorber  is  unpainted  aluminum 
folded  in  deep-V,  accordion  pattern  along  the  vertical  run  of  the  collector. 
The  deep  V is  intended  to  absorb  incident  radiation  through  multiple  reflec- 
tions. Air  flows  on  both  sides  of  the  absorber  effecting  a large  heat  trans- 
fer area.  All  collector  panels  are  joined  at  the  bottom  by  a cold  air  dis- 
tribution duct  and  at  the  top  by  a hot  air  collection  duct.  The  panels  are, 
thus,  all  connected  in  parallel  and  flow  is  from  bottom  to  top.  There  is  a 
manual  summer  vent  and  also  a manual  damper  to  connect  the  collectors  to  the 
house  for  natural  convection  heating. 

The  rock  bin  heat  storage  is  built  of  concrete  block  and  is  located  im- 
der  the  mechanical  room  ( see  Figure  4) . The  exterior  of  the  bin  is  insulated 
with  2-inch  rigid  foam;  the  base  is  not  insulated.  The  bin  is  filled  with 
1-|- inch  gravel  and  has  a volume  of  16,7  m^  (600  ff^).  A galvanized  DHW  pre- 
heat tank  is  buried  in  the  rock  bin.  This  tank  was  not  connected  at  the  time 
of  the  monitoring. 

5.0  PERFORMANCE  MONITORING 

5. 1 Flow  Measurements 

Duct  flow  measurements  were  made  in  January  1979  as  the  initial  step  in 
the  monitoring  data.  At  least  five  holes  were  drilled  in  the  ducts  at  sev- 
eral test  sections,  A hot-wire  anemometer  probe  was  inserted  in  each  hole 
to  three  predetermined  depths  and  readings  of  the  mass  velocity  were  taken. 

On  the  large  ducts,  a total  of  fifteen  measurements  were  made  on  the  cross- 
section.  These  velocities  were  averaged  and  multiplied  by  the  flow-area  to 
calculate  the  volximetric  flow  rate. 
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Typical  results  are  shovm  in  Figures  7 and  8.  These  drawings  show  the 
locations  of  the  measurement  points  and  the  velocity  measured  at  the  point. 
Note  that  the  ducts  have  insulation  on  the  inside  which  is  2.5  cm  (l  in) 
thick.  Most  of  the  test  sections  were  located  in  a section  of  duct  that 
was  straight  for  about  2 m (6  ft)  and  the  flow  was  relatively  steady  and 
the  velocity  distribution  smooth.  The  accuracy  of  these  measurements  was 
estimated  to  be  1 5^, 

The  flow  measurements  of  Figures  7 and  8 are  summarized  in  Figure  9* 

The  mode  of  the  system  was  "collector  heat  direct  to  house".  The  measure- 
ments show  rather  large  imbalances  in  flow.  The  flow  out  of  the  collectors 
is  about  2-|-  times  the  flow  in  the  inlet  duct.  The  flow  entering  the  plenum 
under  the  house  is  about  three  times  the  flow  returning  in  the  cold  air 
duct.  We  repeated  these  measurements  several  times  with  similar  results, 

Mr,  Mattson  independently  measured  the  flow  rates  in  some  of  the  same 
ducts  with  different  instrumentation  and  his  results  shows  similar  im- 
balances in  the  input  and  output  flows.  He  subsequently  checked  the  summer 
vent,  collectors  and  duct  system  for  major  leaks.  Since  no  obvious  leaks 
were  found,  it  was  assiimed  that  the  leakage  must  be  distributed  among  the 
numerous  joints,  seams  and  connections  in  the  system.  Some  of  these  leaks 
are  apparently  within  the  mechanical  room  as  the  pressure  in  this  room  is 
lowered  when  the  fan  is  on. 

There  is  leakage  past  all  the  dampers  which  complicates  the  real  flow 
distribution.  The  crawl  space  has  some  observable  openings  which  let  in 
light  and  must  let  air  out.  Mr,  Mattson  and  Energy  Alternatives  are  cur- 
rently seeking  solutions  to  these  problems. 

5*2  Data  Analysis  Flow  Approximations 

It  is  impossible  to  accurately  and  conclusively  define  the  performance 
of  this  system  due  to  the  leaks  discussed  in  the  previous  section.  The  col- 
lector clearly  is  absorbing  solar  heat,  some  of  which  is  delivered  to  the 
house  and  the  storage.  The  variety  of  immeasurable  leaks  within  the  qystem 
(dampers)  and  between  the  system  and  its  enclosure  and/or  the  ambient  atmos- 
phere obviate  quantitative  description. 

As  a "best  effort"  toward  analyzing  the  data  we  have  made  the  following 
assiunptions: 

(l)  The  average  flov;  rate  of  the  system  is  13.8  m^/min  in  all  modes. 
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(2)  In  the  collector-to-house  mode,  the  output  of  the  collectors  is 
based  on  (1)  and  the  temperature  entering  the  crawl  space  minus 
the  temperature  in  the  collector  inlet  duct. 

(3)  In  the  collector-to-etorage  mode,  the  output  is  calculated  as  in 
(2)  except  the  outlet  temperature  is  taken  as  the  air  temperature 
entering  the  rock  bin  storage, 

(4)  In  the  storage-to-house  mode,  the  temperature  of  the  air  entering 
the  crawl  space  is  taken  as  the  outlet  and  the  cold  air  return  to 
the  rock  bin  as  the  inlet. 

A wood  stove  was  used  to  provide  auxiliary  heat  through  the  monitoring 
period.  A temperature  probe  located  behind  the  wood  stove  reflected  the 
intensity  of  the  fire.  This  temperature  was  used  to  approximate  the  heat 
output  of  the  stove, 

5. 3 Collector  Efficiency  Curves 

Using  the  flow  approximations  described  in  the  previous  section,  the 
efficiency  curves  for  the  Mattson  collector  were  constructed.  This  data  is 
shown  in  Figure  10  along  with  the  HUD  reference  for  this  general  type  of 
collector.  The  Mattson  collector  output  was  defined  as  output  to  the  crawl 
space  and,  thus,  includes  penalties  for  heat  losses  in  the  duct  work.  This 
may  explain  the  apparent  lower  efficiency  in  some  ranges. 

5.4  Hourly  Performance  Summary 

This  section  presents  selected  hourly  and  daily  performance  data  for 
the  Mattson  house  between  April  10  and  Jime  I6,  The  house  was  occupied  for 
most  of  this  time  and  the  solar  system  was  operating.  The  heat  balance  equa- 
tions and  a description  of  the  itemized  data  are  given  in  Figure  11  and  Table 
2.  The  complete  hourly  data  is  given  in  Appendix  I.  A single  table  has  been 
selected  to  familiarize  the  reader  with  the  data  reduction  format. 

Table  3 for  the  day  of  April  30  illustrates  several  system  modes  and 
will  be  discussed  as  an  example.  During  the  early  morning  hours  there  were 
no  large  heat  inputs  and  the  house  temperature  was  slowly  decreasing.  (There 
was  stored  heat  available  from  the  solar  wall  because  the  wood  stove  had  been 
used  on  the  evening  of  the  previous  day.) 

At  1 AM,  the  heat  storage  outlet,  STORE  OUT,  shows  35°C  (95°F)  and  the 
collector  outlet  duct,  COLL  OUT,  shows  26°C  (79°F),  This  indicates  a natural 
convection  back  flow  from  the  storage  bin  to  the  collector  during  the  early 
morning. 
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Between  6 AM  and  11  AM,  the  system  delivers  heat  from  storage  to  the 
house,  STORE  HOUSE,  and  the  house  temperature  increases.  Note  that  the  col- 
lector outlet  temperature  decreases  at  this  time  indicating  leakage  past 
damper  D1. 

At  11  AM,  the  solar  insolation  has  reached  a sufficient  level  and  the 
collectors  start  delivering  solar  heat  to  the  house.  The  solar  radiation  is 
a maximum  at  1 PM  (13  hours)  and  the  collector  outlet  is  a respectable  73°C 
(i63°P)  and  the  air  entering  the  crawl  space  is  63°C  (145°F).  The  average 
air  temperature  in  the  crawl  space  is  30^0  (86°P)  which  would  be  the  air 
temperature  entering  the  house. 

Between  5 and  6 PM,  the  collectors  deliver  a small  amount  of  heat  to 
the  rock  bin  storage.  Between  9 aJid  10  PM,  stored  heat  is  again  delivered 
to  the  house.  The  passive  solar  wall  delivers  heat  to  the  house  d\xring  the 
first  half  of  the  day  and  absorbs  heat  during  the  last  half. 

The  total  heat  input  for  this  24-hour  period  agrees  very  well  with  the 
calculated  heat  load.  This  close  agreement  is  fortuitous  due  to  the  approxi- 
mations involved  in  the  calculations.  The  longer-term  agreement  of  the  aver- 
age of  these  two  quantities  indicate  that  the  many  assumptions  in  the  data 
analysis  are  not  too  faJ*  off. 

Further  scanning  of  hourly  data  in  Appendix  I show  instances  when  the 
control  system  is  not  operating  as  it  should: 

(a)  Tables  AI-2  to  AI-4,  SOLAR  HOUSE  AND  STORE  HOUSE  are  zero  even 
though  house  temperature  is  low* 

(b)  Tables  AI-5  to  AI-7  and  AI-10  to  AI-12.  SOLAR  HOUSE  is  negative 
because  collector  fan  is  running  at  night, 

(c)  Tables  AI-9  to  AI-12.  SOLAR  STORE  shows  that  system  is  running  at 
night  in  this  mode,  probably  cycling  briefly, 

(d)  Tables  AI-22,  AI-26,  AI-29.  System  did  not  collect  heat  even  thou^ 
sun  was  out  and  storage  was  cold, 

6.0  OVERALL  PERFORMANCE  SUMMARY 

The  performance  from  April  10  to  May  12  is  summarized  in  Tables  4 and  5» 
Table  4 is  a siimmary  of  all  data  contained  in  Appendix  I.  These  days  were 
selected  as  best  representing  the  system  performance  as  most  of  the  system  was 
operating  most  of  the  time  and  the  house  was  occupied  most  of  the  time.  Due 
to  the  large  ntunber  of  operational  problems  during  the  monitoring  period,  this 
data  cannot  be  expected  to  represent  the  potential  performance  of  this  system. 
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Table  5 further  condenses  the  data  into  an  overall  summary.  This  data 
shows  that  the  solar  system  delivered  15?^  of  the  incident  radiation  to  the 
house.  Of  the  heat  put  into  the  heat  storage  bin  only  20^  eventually  was 
delivered  to  the  house.  Since  about  half  of  the  collected  heat  went  into 
the  storage,  these  losses  reduce  the  overall  solar  fraction. 

The  second  part  of  Table  5 compares  the  relative  contributions  of  each 
of  the  heat  inputs  to  the  total  heat  required.  The  calculated  heat  load  is 
seen  to  be  about  25^  greater  than  the  measured  total  heat  input.  This 
could  be  due  to  the  conservation  bias  of  the  calculations,  passive  solar 
gain  and/or  assumptions  in  the  data  analysis.  Using  the  total  heat  input 
as  1009^,  the  solar  contribution  is  The  electric  power  and  wood  stove 

are  21^  and  37^  respectively.  During  the  svimmaiy  period  the  average  house 
temperature  was  14«9°C  ( 58°P)  while  the  average  ambient  temperature  was 
5.6°C  (42°F). 

7.0  PROJECTED  ANNUAL  PERF0RI4MCE 

In  this  section  the  data  from  the  monitoring  period  is  used  to  project 
the  annual  performance  of  the  solar  collectors.  Table  6 compares  solar 
radiation  measured  at  the  Mattson  site  with  simultaneous  measurements  at 
the  Bozeman  High  School.  The  readings  agree  to  within  5 "to  10^  for  weekly 
averages.  These  differences  could  be  due  to  atmosphere  variations,  re- 
flections and  shading  variations  and  instrumentation  calibration  tolerances. 

An  overall  collector  performance  curve  was  drawn  from  data  taken  during 
the  monitoring  period.  This  curve  is  used  as  a simple  model  for  converting 
solar  insolation  to  collector  output.  Figure  12. 

Table  7 summarizes  the  estimate  of  annual  performance.  Solar  insola- 
tion data  is  taken  for  the  period  of  July  1978  to  June  1979  a't  the  Bozeman 
High  School.  Degree  day/month  data  is  taken  from  long-term  averages  for 
Bozeman.  Based  on  this  data  and  the  simple  performance  model  of  Figure  12, 
solan  heat  collected  and  heat  loads  are  shown  in  Table  7»  The  annual  solar 
fraction  is  seen  to  be  30^»  Excess  heat  is  available  during  the  summer 
which  may  be  used  for  heating  DHW. 

8.0  GRAPHICAL  PRESENTATION  OF  DATA 

Figures  13  and  14  show  graphs  of  some  of  the  data  for  ten  days  early 
in  May.  The  corresponding  data  tables  for  this  period  are  shown  in  Appendix 
I,  Tables  AI-14  through  AI-24. 
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The  house  air  temperature  is  seen  to  fluctuate  in  between  I5  and  20*^C 
( 59  and  68^P) . The  passive  wall  temperature  follovjs  the  general  trends  in 
house  temperature  but  is  much  smoother  due  to  the  higher  thermal  mass.  It 
appears  that  the  passive  wall  exchanges  heat  primarily  with  the  house  air. 

At  this  time  of  year  the  wall  receives  little  solar  radiation  through  the 
clerestory  windows. 

The  graph  of  crawl  space  temperature  shows  peaks  when  the  solar  sys- 
tem or  the  auxiliary  furnace  delivers  hot  air  to  this  space.  These  peaks 
correspond  to  the  peaks  of  solar  insolation  during  the  day. 

Note  that  on  May  9?  Figure  14f  the  solar  system  did  not  operate  even 
though  there  was  adequate  solar  radiation.  This  may  have  been  due  to  a 
control  malfunction.  The  wood  stove  vias  on  for  most  of  this  day  (see  Appen- 
dix I,  Table  AI-22)  and  the  passive  wall  temperature  increased.  The  crawl 
space  temperature  during  this  period  shows  a small,  general  decrease  in 
temperature.  A slight  rise  in  crawl  space  temperature  corresponds  to  the 
peak  in  house  air  temperature  during  late  afternoon. 

There  is  no  insulation  in  the  floor  over  the  crawl  space.  The  graphs 
show  that  there  will  be  heat  lost  from  the  house  into  the  crawl  space  at 
times  when  the  furnace  or  the  solar  system  is  not  delivering  hot  air  to  the 
space.  The  equilibrium  temperature  of  the  crawl  space  is  around  5°C  (9^P) 
belov;  house  air  temperature. 

The  owner  has  changed  the  design  of  the  heating  system  since  this  per- 
formance monitoring  period.  Ducts  have  been  installed  beneath  the  house  and 
the  crawl  space  is  no  longer  used  as  a plenum.  The  owner  thinks  the  duct 
system  is  an  improvement. 

9.0  CONCLUSIONS 

The  Mattson  solar  system  is  still  in  the  initial  or  shake-down  phase 
of  operation.  The  scope  of  the  data  base  from  this  monitoring  effort  is 
correspondingly  limited.  Major  system  problems  enco\intered  were  (a)  large 
leakcLge  flows,  (b)  control  malfunctions,  and  (c)  a bumed-out  fan  motor. 

Mr.  Mattson  has  solved  some  of  these  problems  and  is  currently  taking  steps 
to  improve  the  system.  The  final  performance  of  the  modified  system  should 
be  much  better  than  the  measurements  and  data  contained  in  this  report. 
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TABLE  1 

MATTSON  HOUSE  HEAT  LOAD 


R 
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U X A 

(Btu/hr  ft^  °P) 

(sq.  ft.) 
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2287 
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Floor 
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.05 

2092 

104.6 

W indows 

1.5 

. 66 
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Skyli^ts 
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*Infiltration: 
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Total  Heat  Load 
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(13)  Temperature  of  air  entering  crawl  space  from  system  (average  of  3 probes) 

(14)  Average  temperature  of  crawl  space  air  (3  probes). 

(15)  Ambient  air  temperature. 
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TABLE  4 

DAILY  PERFORMANCE  SUMMARY,  MATTSON  HOUSE 


Da.y 

SolcLT 

Inso. 

Solar 

House 

Solar 

Store 

Store 

House 

Elect 

Pov;er 

Stove 

Power 

Total 

Input 

Heat 

Load 

House 

Temp 

Ambnt 

Temp 

4/10 

163 

0 

12 

0 

12 

20 

42 

280 

10.1 

- .1 

4/11 

899 

0 

260 

0 

34 

106 

137 

331 

10.9 

.8 

4/12 

1040 

0 

288 

0 

35 

0 

36 

323 

10.8 

1.0 

4/13 

359 

0 

45 

0 

19 

146 

166 

291 

11.4 

2.6 

4/14 

920 

157 

65 

0 

40 

0 

196 

242 

11.9 

4.7 

4/15 

555 

91 

11 

1 

57 

0 

147 

213 

12.9 

6.7 

4/16 

530 

86 

28 

0 

66 

0 

152 

139 

13.3 

9.4 

4/17 

973 

0 

265 

0 

32 

0 

30 

104 

15.1 

12.3 

4/18 

303 

0 

49 

0 

19 

0 

23 

313 

13.0 

4.0 

4/19 

836 

18 

145 

1 

42 

0 

62 

366 

12.2 

1.3 

4/21 

1012 

242 

10 

1 

56 

0 

294 

302 

13.6 

4.4 

4/22 

683 

93 

49 

0 

36 

179 

307 

239 

13.8 

5.3 

4/30 

897 

207 

9 

68 

79 

0 

353 

348 

19.3 

9.0 

5/1 

379 

0 

76 

32 

54 

99 

188 

419 

17.9 

5.5 

5/2 

1080, 

232 

0 

53 

85 

262 

629 

481 

19.0 

4.6 

5/3 

1057 

282 

3 

1 

75 

0 

357 

371 

19.2 

8.1 

5/4 

267 

67 

4 

2 

70 

25 

166 

237 

17.5 

10.6 

5/5 

324 

2 

2 

2 

57 

134 

194 

219 

18.7 

12.4 

5/6 

316 

0 

0 

0 

51 

297 

348 

446 

18.0 

5.0 

5/7 

922 

226 

3 

52 

60 

0 

341 

410 

16.8 

4.5 

5/8 

504 

68 

1 

39 

56 

451 

613 

520 

16.7 

1.1 

5/9 

928 

0 

0 

0 

39 

458 

494 

528 

18.5 

2.5 

5/10 

744 

181 

2 

1 

73 

84 

341 

417 

17.8 

5.2 

5/11 

658 

155 

4 

1 

69 

0 

227 

311 

17.1 

7.9 

5/12 

676 

165 

1 

5 

68 

0 

237 

255 

17.7 

10.3 

TOTALS 

17025 

2272 

1332 

259 

1284 

2261 

6080 

AVERAGE 

8105 

TEMP 

14.9 

5.6 

12 


TABLE  5 

OVERALL  SUMMARY  FOR  25  DAYS* 


MegaJoules 

10^  BTU 

$ 

(1)  Total  solar  insolation 

17,025 

16.1 

100 

(2)  Solar  directly  to  house 

2,272 

2.15 

13 

(3)  Solar  to  storage 

1,332 

1.26 

8 

(4)  Recovered  from  storage 

259 

0.25 

1, 

(5)  Overall  solar  efficiency  (2)+(4)/(l) 

15 

(6)  Storage  efficiency  (4)/  (3) 

20 

Calculated  house  heat  load 

8,105 

7.68 

Total  heat  input 

6,080 

5.76 

100 

Solar  input  [(2)+(4)] 

2,531 

2.4 

42 

Electric  povfer  input 

1,284 

1.22 

21 

Wood  stove  input 

2,265 

2.15 

37 

Average  house  temperature 

14.9°C 

58.1°P 

Average  ambient  temperature 

5.6°C 

42®P 

*This  period,  4/1O  to  5/l2,  was  selected  as  being  most  representative  of  the 
operation  of  the  solar  system. 


TABLE  6 

SOLAR  INSOLATION  COIVIPARISON 


Week  Starting 

Bozeman  High  School 
kWhm  ^ 

Mattson 

kWhm“^ 

3.7 

3.9 

4.0 

4.1 

4/29/79 

4.0 

3.8 

5/6/79 

(3.6) 

3.9 

PROJECTED  YEARLY  PERPORI-IANCE  SUMMARY 
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Figure  1:  Mattson  Solar  House 
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Figure  2;  Mechanical  Equipment  Room 
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Figure  3*  Floor  Plan  Mattson  House 
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Figure  4:  Section  - Mattson  House 
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Figure  5*  Shading  Diagram  Mattson  House 
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Figure  6:  System  schematic  with  transducer  layout  and  control  matrix 
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Figure  8;  Flow  to  and  from  house 
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Figure  Schematic  with  flow  rates  in  cubic  meters/minute 
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Figure  10:  Collector  Efficiency  Curves 
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3000RE«  3UBR  +*+  COrtPUT!:  HOURLY  SUiMS  + + ^ 

3005  '"--  SUMMING" 

3010  N0^N0+1\R£M  HHRS  ' - 

3012  REM  ALL  THE  EOLLOUING  ARE  IN  MEGA  JOULES 

3015  S(1  ) = (li(1  ) + . 041*48. 12  + 3.6  + . 95  \REM  TOTAL  SOLAR  INPUT  CORRECTEIi 
3020  S(2)=D(3)*D(4)*13.8  + 1 .01 +1  .2*.06KD{16)-I!(15)  )\REM  SOLAR  TO  HOUSE 
3025  S(3)  = (1-D(4))+D(3)t13. 3*1. 01  + 1. 2+.06*(D(17)-j[)(15)?\R£M  SOLAR  TO  STORAGE 
3030  S(4)-(1-!T(3))+ri(45  + 13.8*1 .01  + 1 .2+.06  + 1D(16l-[i(21  ) )\REM  STORAGE  TO  HOUSE 
3035  S(5)=D(2)*3.6*.9  \ REM  ELECTRICAL  POL'ER , CORRECTEIi 

3040  S(6)  = 1 .2+{D(12)-[i(9>)\  IF  S(6)<3  THEN3041  ELSE3045\  REM  STOVE  OUTPUT 

3041  S(6)=0\REM  ZERO  STOVE  FOR  SMALL  AND  NEGATIVE  VALUES 
3045  S(7)=1 .8+(T1-D(22))\REM  PASSIVE  STORAGE  UALL  CHANGE 
3050  S(8)=S(2)+S(4)+S(5)+S(6)+S(7)  \R£M  TOTAL  INPUT  ENERGY 
3055  S(9)=1.37+(D(9)-D(11 ))\R£M  HOUSE  HEAT  LOAD 

3059  REM  PUT  IN  TEMPERATURES 

3060  S ( 1 0 ) - D ( 9 ) \ S (1 1 ) -Li  ( 1 9 ) \S  ( 1 2 ) -D  ( 1 7 ) \S  n 3 ) -^D  ( 1 6 ) \S  ( 1 4 ) -=0  ( 1 8 ) \S  ( 1 5 ) =0  ( 1 1 ) 
3065  T1=D(22)\REM  SAVE  STORAGE  UALL  TEMPERATURE  FOR  DIFFERENCE 

3130  !#1  ,23I,R,Z4I,M0,X3I,[i0,3:5I,H0,TAP(15), 

3140  FOR  1=1  TO  15\T(I )=T(I)+S(I)\REM  SUM  FOR  DAY.  PRINT  HR  DATA 
3150  'ttl  .%7F1  ,S(I).\NEXT\*tt1 
3160  R£TURN\'lt1 

3200REM  * + + CALC  AVERAGES  3.  FINAL  QUANTITIES 

3210  FOR  1=10  T015\  REM  AVERAGE  TEMPERATURES  (NOTE  HEATS  NOT  INCL.,) 

3220  T(I)=T(I)/NO\R£M  DAILY  AVGS  = SUM  / NR  HRS 
3230  NEXTXRETURN 

3400REM  SUBR  * + +•  PRINT  REPORT  HEADING 


3405 

!t{1\'ftr'  DAILY  PERFORMANCE 

SUMMARY  FOR 

MATTSON 

HOUSE 

(ID=-" 

,R$,"')" 

3410 

•iilX'KI  ,"REC«  MO  DAY 

HR  " . 

TAB(16) 

« 

3420 

‘!t1,"  SOLAR  SOLAR 

SOLAR 

STORE- 

ELECT 

STOVE 

PASS 

TOTAL 

HEAT  ", 

3430 

'«1,"  HOUSE  COLL 

STORE 

HOUSE 

CRAUL 

AMBNT  " 

3440 

•ffl  ,TAB(16) , 

* 

3450 
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Figure  11:  Computer  Heat  Balance  Ecjuations 


25 


aaioaiioo  ivioj, 


Figure  12:  Overall  Collector  Performance  Ci 
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Figure  13:  Graph  of  5 Days  Data  in 


Figure  14;  Graph  of  5 Days  Data  in  May 
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TABLES  OF  HOURLY  DATA  FOR  MATTSON  HOUSE 
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APPENDIX  II 


PASSIVE  OPERATION  OP  SOLAR  COLLECTORS,  MATTSON  HOUSE 


II -1 


PASSIVE  OPERATION 

During  the  early  part  of  the  monitoring  period,  the  Mattsons  had  not 
moved  into  the  house  and  the  solar  collectors  were  being  operated  in  the 
natural  convection  or  passive  mode.  The  natural  convection  damper,  Figure 
6,  was  open  and  the  fan  was  off.  The  fan  motor  was  burned  out  during  part 
of  the  monitoring  period  and  the  system  was  run  in  the  passive  mode. 

Figure  AII-1  shows  the  results  of  a single  measurement  of  the  collec- 
tors* output  in  the  passive  mode.  This  test  was  done  shortly  after  solar 
noon  on  a veiy  cleax  day  neair  the  end  of  January, 

The  flow  rate  out  of  the  duct  was  5*93  m^/min  (207  cfm)  and  the  outlet 
temperature  was  105*^0  (221^F).  The  incident  solar  radiation  in  the  plane 
of  the  collector  was  49»43  kW  (I5f670  BTU/hr)  and  the  collectors  were  de- 
livering 11,2  Wnf  (3550  BTU/hr)  to  the  house.  The  conversion  efficiency 
was  23^. 

The  computer  heat  flow  equations  were  modified  for  the  passive  mode, 
Figure  AII-2.  The  passive  collector  output  was  calculated  using  a constant 
flow  rate  of  6 m^/min  and  the  inlet  and  outlet  duct  temperatures. 

This  program  produced  the  results  shown  in  Tables  AII-1  through  All -7 
for  the  first  week  of  April.  The  overall  efficiency  of  the  collector  dur- 
ing this  week  was  around  109S.  The  collector  output  temperature  was  a maxi- 
mum of  46. 5°C  (11 6°F)  on  April  8,  much  lower  than  the  January  measurement, 

I suspect  that  some  of  the  ductwork  or  dampers  were  open  in  the  mechanical 
room  which  was  diluting  the  output  temperature.  During  the  early  part  of 
the  monitoring  period,  the  damper  blade  of  D1  was  slipping  on  its  shaft  and 
would  sometimes  be  opened  20°  or  more. 

The  passive  heat  was  sufficient  to  maintain  the  average  house  tempera- 
ture around  10°C  (l8°F)  above  ambient  conditions  during  this  time.  This 
heat  came  entirely  from  the  passive  windows  and  the  collectors. 
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51  cm 
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cm 
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< - v=.  96“^/ sec 


Seimple  Calculations;  '^:00  pm,  Jainuary  ?6,  1979 


Duct  Area  * 0.103m^ 

Duct  Flov;  = 0,96  m/sec  x 0.103  x 60 
= 5*93  m3/min 


INC  I DMT  POWER 


ABSORBED  POWER 


1.02  kW/m?  X 48.12  m2 

49.43  kW 

5.93  ^ 1.024  X 1.2 

60 

11.2  kV/ 


.1  30c ) 


EFFICIMCY  = 11.2  x 100  = 22,efo 

49.4 


FIGURE  AII-I1  PASSIVE  PLOW  PROFILE  AND  CALCULATIONS 


II 


i.^OOCREH  SUBR  *++  COMPUTE  HOURLY  SUMS  :*t«^++t.**++ttft+*++*++t+tt+tt+ttt+t 
3005  SUMHING" 

3010  NO  = NO  + 1\RE,*1  «MRS 

3012  REM  ALL  THE  FOLLOUING  ARE  IN  NEGA  JOULES 

3015  S(1  ) = (!!(  1 ) + .04)  MS.  12t3.A+.95\REM  TOTAL  SOLAR  INPUT 

3020  S(2>=  6 + 1 .01+1 .2+.06+{I.i(19)-D<rj)  )\REM  SOLAR  TO  HOUSE 

3021  IF  S(1)>1  THEN  3025  ELSE  3024 

3024  S(2)=0 

3025  S(3)  = (1-D(4)  ) + [i(3)+13. 8 + 1 .01  + 1 .:?+.()6+(I,i(17)--D(15))\RE«  SOLAR  TO  STORAGE 

3030  S(4)  = n-ri(3))+H(4)  + 13. 8+1. 01*1. 2+. 06+(n(16)-Ii(21  D'nREM  STORjIGE  TO  HOUSE 
3035  S(5)=D(2)+3.6+.9  \ REN  ELECTRICAL  POUER 

3040  S(6)  = 1 .2+(ri(12)-D(9)  )\  IF  S(6)<3  THEN3041  ELSE3()45\  REN  STOUE  OUTPUT 

3041  S(6)=0\REM  ZERO  STOVE  FUR  SHALL  AND  NEGATIVE  VALUES 
3045  S(7)  = 1 .8+(T1-[!(22)  )\REN  PASSIVE  STORAGE  UALL  CHANGE 
3050  S(8)rS(2)+S(4)+S(5)+S(6)+S(7)  \REN  TOTAL  INPUT  ENERGY 
3055  S(9)-1 .37*(D(9)-D( 1 1 ) )\REN  HOUSE  HEAT  LOAD 


3059  REN  PUT  IN  TEMPERATURES 

3060  S<10)=B(9)\S(1 1 ) = B(19)\S(12)^ri(17)\S(13)-:D(16)\Sn4)-D<18)\S(10)=U(1 1 ) 
3065  T1=D(22)\REM  SAVE  STORAGE  UALL  TEMPERATURE  FOR  DIFFERENCE 

3130  Ml  ,23I.R,%4I.M0,X3I,D0,J:5I,H0,TAB(15>. 

3140  FOR  1^1  TO  15\T(I)=T(I)+S(I)\R£M  SUM  FOR  DAY,  PRINT  HR  DATA 
3150  Ml  ,X7F1  ,S(I),\NEXT\'t}1 
3160  RETURN\M1 

3200REM  *++  CALC  AVERAGES  Jt  FINAL  OUANTITIES 

3210  FOR  1^-10  T015\  REM  AVERAGE  TEMPERATURES  (NOTE  HEATS  NOT  INCL...) 

3220  T(I)=T(I)/N0\R£M  DAILY  AVGS  - SUM  / NR  HRS 
3230  NEXTXRETURN 

3400REM  SUBR  +++  PRINT  REPORT  HEADING 


3405 

MIXMI"  DAI 

LY  PERFORMANCE 

SUMMARY  FOR 

MATTSON 

HOUSE 

(TD=-" 

vR« 

3410 

M1\M1  ,"RECH 

MO  DAY 

HR  ", 

TAB(16) 
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Ml,"  SOLAR 

SOLAR 

SOLAR 

STORE 

ELECT 
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TOTAL 

HEAT 
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CRAUL 

AMBNT  " 

3440 

M1,TAB(16), 

3450 

Ml,"  INSO 

HOUSE 

STORE 

HOUSE 
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POUER 

UALL 

INPUT 

LOAD 

3460 

Ml.  " TEMP 
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OUT 

IN 

SPACE 

TEMP" 

3464 

Ml  ,TAB(16) , 

3465 

Ml,"  MJ 

MJ 

MJ 

M.J 

MJ 

M j 

MJ 

MJ 

MJ 

3466 

Ml,  " C 

C 

C 

C 

C 

C" 

3495 

RETURN 

3500REM  SUi^R  * + + PRINT  SUMMARY  DATA 


3510  M1\M1  ."DAY".X4I.M1  ^:3T,D1  , :i5I  ..NO  ..TAB(  1 5) , 
3520  FOR  1=1  TO  15 


3530  M1,%7F1,T(I).\.NEXT\M1,"  PASSIVE  "\M1 
3535  FOR  I = 1T0  10\M1  \ NEXT 
3540  RETURN 
3550END 


FIGURE  AII-2|  PASSIVE  EQUATIONS 


DAILY  PERFORfJANCE  SUMMARY  FOR  MATTSOR  HOUSE  (ID= 
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ABSTRACT 

This  paper  summarizea  the  thermal  perform- 
ance of  a passive  solar  house.  Hourly  data 
taken  over  a three  month  period  during  one 
of  the  coldest  winters  on  record  showed 
that  the  solar  system  was  providing  65^  of 
the  heating  load.  The  average  house  temp- 
erature was  19°  C and  the  average  daily 
variation  was  6.4°  C.  A low-cost  data  ac- 
quisition system  utilizing  a Radio  Shack 
TliS-80  computer  was  developed  for  this 
project. 


1.  INTRODUCTION 

We  originally  intended  to  present  perform- 
ance data  on  a passive  solar  house  and  a 
nearby  active  air  solar  house  in  this  paper. 
During  the  test  period  the  recently  com- 
pleted active  house  experienced  a fan  fail- 
ure as  well  as  some  control  problems.  Plow 
measurements  indicated  very  large  imbalemces 
of  flow  in  both  the  collector  and  house  loop 
apparently  due  to  air  leaks  in  the  collect- 
ors and  dampers.  Because  of  these  problems 
the  active  house  data  was  judged  to  be  of 
limited  value  and  is  not  presented  in  this 
paper. 


2.  i'ASSIVi;  HOUJE  DL^SCHIFTIOll 

The  passive  solar  iiouse  is  located  in 
Bozeman,  MT  at  cui  elevation  of  14b0  m in  a 
broad  mountain  valley.  Environmental  data 
are  shown  in  Table  1.  The  house  was  design- 
ed and  built  by  the  author  during  1976-77 
and  it  has  been  occupied  by  the  author  and 
his  family  since  September  1977*  A photo- 
graph of  the  house  is  shov/n  in  Figure  1. 

The  house  is  a two-stoiy  v;ood-framed  struc- 
ture with  a usiible  lloor  area  of  about  19O 
m^.  The  walls  are  framed  with  5^15  c"  studs 
on  60  cm  centers  and  the  roof  is  spanned 
by  5^25  cm  purlins.  The  walls  and  the  floor 
are  insulated  with  15  cm-tliick  glass  fiber 
batts  and  the  roof  with  25  cm-thick  batts. 

A continuous  vayioi  barrier  of  polyethylene 
\ja5  used  inside  all  insulated  surfaces. 


Table  I;  Bozeman  Climatological  Data  (I) 


Deg.  Day 

Deg. Day 

Insol. 

Insol. 

Avg. 

1978-79 

1977-78 

, 1978-79 

°C-d 

°C-d 

kWh/m 

(60°) 

Oct 

322 

_ 

4.41 

5.24 

Nov 

542 

- 

1.99 

2.92 

Dec 

682 

864 

1.93 

2.35 

Jan 

762 

1019 

2.53 

2.8 

Feb 

614 

730 

3.81 

3.48 

Mar 

596 

- 

4.98 

3.97 

Pig.  1 : Passive  solar  house 


Fenestration,  excluding  the  collector  glaz- 
ing io  12^  of  the  floor  area  and  all  windows 
are  double  glazed.  Electricity  is  used  for 
lighting,  hot  water  and  all  appliances. 

2. 1 Passive  Colar  System  Description 

Radiation  enters  a vertical,  double-glazed 
aperture  viluch  faces  due  south.  The  total 
area  of  the  apertuie  is  43.6  m^.  The  center 
4m^  is  not  backed  by  the  solar  wall  and 
serves  as  a v/indow. 

The  solar  wall  consists  of  20  vertical  tubes 
filled  with  water.  The  tubes  are  filament 
v;oujid  FRP,  40  cm  in  diameter,  5*2  cm  long 
and  contain  674  kg  of  water  each.  The  total 
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Btorage  mass  of  the  wall  is  about  13|430  kg. 
The  tubes  are  painted  flat  green  on  the 
portion  facing  the  sun  and  are  spaced  3«3  on 
apart.  Heat  is  transferred  from  the  solar 
wall  to  the  house  by  natural  convection  and 
radiation. 


The  calculated  average  heat  load  for  the 
house  using  ^ air-change/hour  is  1.08  Mj/h°C 
(1.41  Mj/h®  C with  insulating  curtain  open; 

0.9  Mj/h”  C closed).  The  heat  load  for  the 
house  with  the  collector  wall  replaced  by  a 
conventional  wall  is  O.9O  Mj/h°  C. 


Design  calculations  indicated  that  a movable 
insulation  system  would  be  needed  to  control 
the  heat  loss  through  the  collector  glazing- 
at  night.  The  insulation  system  selected 
for  this  house  consists  of  two  thick  fabric 
curtains  which  resemble  window  shades, 
Figures  2 and  3.  The  curtains  are  mounted 
on  a common  shaft  and  actuated  by  a geared 
electric  motor.  The  outside  of  the  curtain 
was  fabricated  from  6 oz.  white,  dacron 
sailcloth  aind  the  inside  from  6 oz.  nylon 
"pack  cloth"  having  a urethane  coating.  The 
insulation  consists  of  4 layers  of  10  oz. 
Polarguard,  a continuous  polyester  fiber. 

The  thickness  of  the  curtain  is  about  12  cm. 
Sewn  along  both  vertical  edges  of  the  cur- 
tain is  a bolt  rope  which  runs  in  a contin- 
uous plastic  guide  attached  to  the  window 
frame  to  minimize  heat  loss  by  convection 
around  the  edges. 


2. 2 Auxiliary  Kood  Furnace 

The  auxiliary  heat  is  supplied  hy  a wood 
burning  furnace.  The  heat  from  the  furnace 
is  normally  directed  to  the  solar  storage 
wall.  The  furnace  always  operates  at  full 
output  with  the  storage  wall  acting  to  mod- 
erate the  amount  of  furnace  heat  entering  the 
house  itself. 

The  furnace  consists  of  a fire  box  and  a 
heat  exchanger.  It  was  designed  by  the 
author  and  fabricated  from  11  and  13  gage 
steel  sheet.  The  fire  box  is  lined  with  fire 
brick  and  is  a base  burning  design.  Pre- 
heated air  enters  the  combustion  chamber 
through  three  adjustable  vents.  The  combus- 
tion gases  leave  the  fire  box  below  the 
grate  and  pass  through  a heat  exchanger.  A 
draft  induction  fan  is  used  when  starting  a 
new  fire.  Heat  is  moved  from  the  furnace  to 
the  solar  wall  with  an  air  loop  powered  by 
a conventional  furnace  blower  and  thermostat. 
The  only  hot  air  outlets  to  the  house  are  at 
the  base  of  the  solar  storage  wadi.  The  cold 
air  returns  from  a manifold  located  on  the 
other  side  of  the  solar  wall. 

If  the  house  temperature  drops  to  around 
15®  C and  cold  and  cloudy  weather  is  pre- 
dicted a fire  is  built.  The  fire  box  is 
filled  with  about  75  of  pine  and  ignited. 
The  furnace  reaches  full  output  during  the 
next  half  hour  and  bums  for  about  6 hours. 

A second  load  of  fuel  is  sometimes  added  to 
the  fire  box  in  extreme  weather. 


Fig. 2:  Storage  wall  and  curtain 


1.  Water  storaige  wall 

2.  Glazed  aperture 

3.  Insulating  curtain 

4.  Auxiliary  hot  air 

5.  Cold  air  return 

6.  First  floor 

7.  Second  floor 


Pig.  3:  Schematic  of  passive  solar  system 


3.  TEST  RESULTS 
3. 1 Introduction 

The  data  acquisition  system,  transducers  and 
data  accuracy  are  discussed  in  the  Appendix. 
Hourly  data  was  processed  to  calculate 
hourly  average  powers  and  temperatures.  The 
basic  balance  equation  was  Power  Input  = 
Power  Output  + Storage  Power.  Four  Power 
inputs  were  considered:  solar  power  on 
storage  wall,  solar  power  directly  into 
house,  furnace  power  into  house  and  electri- 
cal power  dissipation  in  the  house. 

Solar  insolation  was  measured  on  a vertical 
plane  inside  the  glazing.  Solar  power  on  the 
storage  wall  is  calculated  by  multiplying 
this  Insolation  by  the  aperture  area  and  the 
fraction  of  time  the  insulating  curtain  was 
open  during  that  hour.  Solar  power  directly 
into  the  house  wan  calculated  from  the  inso- 
lation times  the  area  of  the  south  facing 
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windows.  Purnaoe  power  was  caloulated  from 
the  inlet  and  outlet  temperatures  and  flow 
rate.  Electrical  dissipation  was  estimated 
as  0.6  of  the  total  electrical  power  coming 
into  the  house. 

Storage  power  was  considered  from  two 
sources:  energy  changes  in  the  solar  wall 
and  energy  changes  in  the  house  itself. 

Solar  wall  power  was  calculated  from  the 
mass,  specific  heat  and  average  temperature 
change  of  the  water  per  hour}  52  Mj/°  C x 
A T/h.  The  house  storage  power  was  estima- 
ted as  8 Mj/°  C X A T/h. 

The  output  power  was  calculated  from  the  six 
quantities  above  and  divided  by  the  tempera- 
ture difference  between  the  house  and  ambient 
air  to  produce  a heat  load  factor. 

3.2  Daily  Performance 

Sample  hourly  data  in  Tables  II  and  III 
illustrate  the  thermal  behavior  of  the  sys- 
tem on  two  different  days.  Typically  during 
the  early  morning  heurs  the  temperatures  of 
the  tubes  decrease  as  they  deliver  power  to 
the  house  and  the  house  temperature  reaches 
a minimum  around  sunrise.  The  sun  charges 
the  tubes  and  heats  the  house  during  the 
day  and  both  tubes  and  house  reach  a maximum 
near  sunset.  The  tubes  are  normally  3 to  9°C 
above  the  house.  This  difference  depends  on 
the  heat  load  since  the  heat  is  removed  from 
the  tubes  by  natural  convection.  Similarly 
the  temperature  variation  in  the  house  is 
usually  greater  for  low  ambient  temperatures. 
The  overall  heat  transfer  coefficient  be- 
tween the  tubes  and  the  house  was  around 
21  Kj/hm^  ®C. 

Table  II:  January  8,  very  cold  and  clear 


Table  III:  March  1,  a clear  moderate  day 


TBIE 

SOLAR 

TUBE 

TUBE 

HOUSE 

AMB. 

POWER 

POWER 

TEMP 

TEMP 

TEMP 

h 

Mj/h 

Mj/h 

°C 

°C 

°C 

1 2 am 

0 

-18.4 

25.0 

16.7 

-27.8 

4 am 

0 

-12.7 

23.7 

14.9 

-27.2 

8 am 

0 

-19.9 

22.1 

13.9 

-26.1 

1 2 pm 

109.2 

47.3 

23.0 

18.7 

-16.4 

4 pm 

38.1 

45.1 

28.0 

22.4 

-16.4 

8 pm 

0 

-15.9 

27.0 

19.3 

-24.1 

1 2 am 

0 

-17.5 

25.6 

17.5 

-23.7 

24  Hour 

Total 

24  Hour  Average 

637.1 

32.6 

24.8 

17.7 

-22.5 

In  Table  II  the  24  hour  net  tube  energy  for 
January  8 was  + 32.b  MJ ; enough  reserve 
energy  to  heat  the  house  for  about  two  hours. 
The  previous  6 days  were  similarly  cold  but 
clear  and  the  house  was  ^0(y/}  solar.  During 
this  period  the  daily  collector  efficiency 
was  around  January  9 was  cloudy  and 
the  furnace  was  fired  at  1 1 am  when  the 
house  temperature  dropped  to  15«8°  C. 


TIME 

SOLAR 

TUBE 

TUBE 

HOUSE 

AMB. 

POWER 

POWER 

TEMP 

TEMP 

TEMP 

h 

Mj/h 

Mj/h 

°C 

°C 

°C 

12  am 

0 

- 8.3 

22.4 

17.2 

- 7.7 

4 am 

0 

- 7.6 

21.5 

16.0 

- 7.4 

8 am 

0 

- 5.7 

20.8 

16.3 

- 4.7 

12  pm 

116.8 

52.9 

22.3 

21.0 

- 1.9 

4 pm 

21.6 

24.4 

27.3 

23.5 

- 5.2 

8 pm 

0 

-10.1 

26.7 

20.7 

-10.3 

12  am 

0 

- 9.5 

26.1 

19.6 

-14.9 

24  Hour  Total  24  Hour  Average 

540.4  194.0  23.7  19.2  - 7.0 


From  Table  III  the  24  hour  net  energy  for 
March  1 was  194  MJ  or  enough  to  supply  the 
house  for  another  24  hours  without  sun  or 
auxiliary  power.  The  average  collector  effi- 
ciency was  54^  on  this  day. 

3. 3 Performance  During  Monitoring  Period 

Hourly  data  was  recorded  between  November  28 
to  March  6.  The  heating  degree  days  during 
the  test  period  were  17^  above  average 
creating  a challenging  test  environment.  The 
average  house  temperature  was  19.0°  C with 
an  average  differential  of  3.7°  0 between 
the  upper  and  lower  floor.  The  average  daily 
variation  in  house  temperature  was  6.4°  C. 

During  the  test  period  electrical  power 
dissipated  in  the  house  was  7.7  GJ.  Thirteen 
fires  were  built  in  the  axixiliary  v/ood  fur- 
nace which  put  a total  of  10.3  GJ  into  the 
house  (about  0.8  cord  of  pine  was  consvuned). 
The  total  heat  requirement  of  the  house  was 
52  GJ  based  on  the  experimental  heat  loss 
factor  of  0.83  Mj/h-°C.  Based  on  these 
measurements  the  sun  provided  34  GJ  or  69j5 
of  the  total  heating  energy  during  December, 
January,  and  February. 

Measurements  of  the  temperature  drop  across 
the  insulating  curtain  implied  that  its 
overall  heat  transfer  coefficient  was  4 Kj/ 
hm^  °C.  This  value  is  much  larger  than  anti- 
cipated and  may  be  due  to  convection  caused 
by  inadequate  sealing  around  the  edges  of 
the  curtain. 

The  average  efficiency  of  the  furnace  was 
determined  by  weighing  the  total  wood  put  in 
during  each  firing  and  then  measuring  inlet 
and  outlet  duct  temperatures  and  outlet  flov; 
rate.  The  heating  value  of  the  wood  was 
taiken  as  18.4  Mj/kg.  Using  these  data  the 
average  overall  efficiency  of  the  wood  fur- 
nace was  58.7/  during  11  tost  firings. 

3.4  Projected  Annual  Performance 

During  the  entire  heating  season  a total  of 
16  fires  were  built  in  the  auxiliary  wood 
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furnace,  the  first  on  November  11  and  the 
last  on  January  27.  A total  of  11 40  kg 
(about  one  cord)  of  wood  was  burned  produc-. 
ing  an  estimated  12.7  GJ  of  heat.  The  annual 
electric  power,  dissipated  as  heat  into  the 
house,  contributed  about  15  GJ  (700  kWh/mo). 
Using  an  average  degree-day's  value  of  4484 
the  annual  heat  load  of  the  house  would  have 
been  90  GJ  and  the  annual  solar  fraction, 
6?^.  Since  the  test  period  this  winter  was 
275®  colder  than  noinnal,  the  average  annual 
solar  fraction  is  probably  about  79^. 

The  insulating  curtain  is  kept  in  the 
raised  position  from  May  to  October.  No 
overheating  occurs  during  the  summer  but 
on  warm  and  sxmny  fall  days  excess  heat 
must  be  vented  for  comfort. 


4.  SUMMARY 

The  test  results  show  that  this  system  pro- 
vides a high  percentage  of  the  house  heat 
load  by  solar  conversion,  even  during 
extremely  cold  weather.  It  seems  unlikely 
that  the  overall  solar  fraction  could  be 
significantly  improved  by  economically 
aidmissable  changes  in  either  aperture  area 
or  storage  thermal  capacity  (2,3,4).  The 
conversion  efficiency  of  the  wood  furnace 
and  the  heat  loss  factor  of  the  house  it- 
self could  not  be  greatly  improved.  Reduc- 
tions of  over  105®  in  the  overall  heat  loss 
would  result  if  the  movable  insulation 
could  be  constructed  so  as  to  have  a heat 
transfer  coefficient  equivalent  to  a normal 
wall.  A reduotion  in  the  24  hour  tempera- 
ture variations  would  increase  the  comfort 
of  the  house. 

The  TRS-80  data  acquisition  system  performed 
to  our  satisfaction.  It  provided  a large  and 
densely  sampled  data  base  from  which  to  con- 
struct this  performance  summary. 
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6.  APPENDIX;  INSTRUMENTATION 

A. 1 . TRS-^  Data  Acquisition  System 

The  data  acquisition  and  processing  require- 
ments for  solar  performeince  monitoring  aire 
best  satisfied  by  a computer-based  system 
programmable  in  a hi^  level  language.  Due  to 
budget  limitations  we  developed  our  own  sys- 
tem using  the  maiss-produced  Radio  Shack  TRS- 
80  with  4K  memory  and  Level  II  BASIC.  This 
unit  sells  for  about  $700  and  includes  a 
cassette  tape  for  storing  programs  and  data. 

Ke  developed  two  main  accessories  for  this 
unit;  a battery  operated  power  supply  and  an 
analog/digital  interfewse.  Figure  A1. 

The  power  supply  contains  circuitry  to  con- 
tinuously charge  a conventional  12v  lead-acid 
o&r  battery  when  110  vac  power  is  available. 
Additional  circuitry  converts  the  battery  vol- 
tage into  the  voltages  necessary  to  power  the 
computer  and  the  a/d  interface.  A smedl  bat—  * 
tery  is  sufficient  to  power  the  eystm  for 
about  5 hours  during  a utility  outeige. 

The  a/D  interface  incorjjorateB  a front  panel 
with  plugs  for  connecting  the  transducers,  a 
signal  conditioning  board,  an  a/d  interface 
board  and  a real-time  clock  with  month,  day, 
hour,  minutes  and  seconds  (5).  The  system 
accepts  40  transducer  inputs:  1 solar  insola- 
tion, 2 electric  povier,  5 status  switches,  8 
pulsed  inputs  ( flow  meters) , 8 averaging 
temperature  inputs  (3  probes  each  channel) 
and  16  single  probe  temperatures.  The  signal 
conditioning  and  transducer  power  supplies 
are  Integrated  into  the  chassis.  Analog  vol- 
tages between  1 and  4 volte  are  converted 
into  a 12  bit  number  giving  a resolution  of 
about  1 mv. 

The  data  acquisition  is  structured  with  a 
program  written  in  BASIC.  The  inputs  are 
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Bcarmed,  the  voltages  converted  to  the 
desired  parameter  using  calibration  con- 
stants and  equations  in  the  BASIC  prOigram  and 
the  current  data  dl'aployod  on  the  video 
screen.  Averaging  and/or  other  data  process- 
ing is  then  executed  and  the  clock  is  read. 
At  the  end  of  the  averaging  period,  speci- 
fied in  BASIC,  the  data  is  v/ritten  onto 
cassette  tape. 


I'ig.  A1;  TP.S-80  Data  Acquisition  System 


The  system  used  for  monitoring  the  passive 
house  completed  about  ^00  scisis  per  liour 
and  wrote  hourly  averages  on  tlic  tape,  liach 
side  of  an  audio  quality  tape  held  about 
11,000  numbers.  Data  on  these  tapes  \jere 
transferred  to  diskette  for  processing  on  a 
laiger  computer.  He  have  been  using  three  of 
these  TOS-80  data  acquisition  systems  for  5 
months  to  monitor  a variety  of  solar  heating 
projects. 


A. 2.  Transducers  and  Accuracy 

Temperatures  were  measured  using  the  AD-590 
made  by  Analog  Devices.  A comouter  pro.'rram 
is  used  with  the  data  system  to  calibrate 
the  probes.  The  absolute  accuracy  of  the 
temperatures  was  better  than  1°  C while  the 
accuracy  of  small  temperature  changes  was  at 


least  0.1°  C. 

The  solar  tremeducer  was  positioned  behind 
the  glaelng  In  the  center  of  the  solar  wall 
to  measiire  actual  solar  input  to  the  storage 
wall.  The  transducer,  based  on  a silicon 
cell,  has  been  referenced  to  the  NOAA  standard 
and  has  an  accuracy  of  5^.  Electric  power 
input  was  monitored  by  measuring  the  current 
in  the  input  mains  with  a clamp-on  ammeter. 
This  transducer  is  cailibrated  on-eite  using  a 
watt  meter  euid  has  an  accuracy  of  at  least 
yfL,  The  status  inputs  used  a switch  or  relay 
to  record  the  status  of  the  fan  and  the  insu- 
lating curtain. 

The  furnace  input  was  determined  from  the 
temperature  difference  in  the  supply  and  the 
return  air  ducts  and  the  flow  rate.  The  accu- 
racy of  this  measurement  is  assumed  to  be  105^. 
Air  leaks  located  beyond  the  measurement  sta- 
tion could  further  increase  the  uncertainty 
of  the  actual  furnace  power  input.  The  weight 
of  wood  burned  in  the  auxiliary  furnace  was 
determined  to  within 

A total  of  27  temperature  probes  were  used. 
Four  sets  of  three  temperature  probes  were 
immersed  in  the  water  in  the  centers  of  the 
heat  storage  tubes.  These  probes  were  placed 
at  varying  levels  in  an  attempt  to  compensate 
for  stratification  and  to  arrive  at  em  average 
storeige  temperature.  The  temperature  change  in 
the  tubes  is  on  the  order  of  1°  C per  hour. 

A temperature  differential  of  5*^  C may  exist 
between  the  top  and  the  bottom  of  a tube  at 
any  given  time.  Analysis  of  the  data  showed 
that  short  term  average  power  calculations 
were  noticeably  effected  by  the  location  and 
number  of  temperature  samples.  Daily  average 
power  calculations,  however,  were  relatively 
independent  of  the  probe  locations. 
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16  copies  of  this  pubiic  document  were  published  at 
an  estimated  cost  of  $9.25  per  copy,  for  a totai  cost  of 
$148.00,  which  includes  $148.00  for  printing  and  $.00 
for  distribution. 


